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STEADY  STATE  NUMERICAL  SOLUTION 
OF  MAGNETICALLY  INSULATED  CHARGE  FLOW  IN 
COAXIAL  GEOMETRY 

I.  INTRODUCTION 

Impressive  results  have  been  achieved  experimentally  with  the  production  of  intense  proton  fluxes 
in  radial  diodes.24  The  high  ion  production  efficiencies  of  these  devices  make  them  attractive  candidates 
for  use  as  driver  sources  in  future  ion  beam  inertial  confinement  fusion  systems.  In  addition,  magneti¬ 
cally  insulated,  coaxial  vacuum  feed  lines5  8  form  an  integral  part  of  pulsed  power  systems  for  a  wide 
variety  of  applications.  Clearly,  a  computer  code  capable  of  predicting  current  flow  in  coaxial  systems 
for  realistic  values  of  operational  parameters  would  be  a  valuable  tool.  At  the  present  time  the  only 
relativistic,  two-specie,  cylindrical  geometry,  theoretical  formulation  using  cycloidal  elec'.on  orbits  is 
that  of  K.  D.  Bergeron.1  Its  foundation  rests  on  the  steady-stale  magnetic  insulation  model  formulated 
independently  by  Sudan  and  Lovelace,9  10  and  at  about  the  same  time  by  Ron,  Mondelli,  and  Ros- 
toker.11  In  this  model,  electrons  form  a  cloud  near  the  cathode  in  which  each  electron  traces  an  indivi¬ 
dual,  single-arc  trajectory  which  begins  and  terminates  on  the  solid  cathode  surface  (See  Fig  1.)  All  of 
these  electrons  see  the  full  steady-state,  self-consistent  electric  and  magnetic  fields  during  their  entire 
orbit  and  all  orbits  are  identical.  There  exists  a  competing  magnetic  insulation  model  formulated  by 
Antonsen  and  Ott12  and  foreshadowed  in  the  work  of  Creedon'1  and  Buneman.14  This  model  assumes 
that  all  electrons  in  the  cathode  sheath  have  been  emitted  from  the  cathode  surface  onto  electrostatic 
equipotential  contours  during  an  adiabatically  rising  diode  voltage  pulse.  In  the  equilibrium  state,  they 
F.  x  B  drift  along  their  respective  contours  in  a  Brillouin15  type  flow,  never  again  intersecting  the 
cathode  (See  Fig.  I).  This  formulation  is  based  on  a  vanishingly  small  Larmor  radius  approximation. 
The  cylindrical  geometry  E  x  Adrift  model  has  recently  been  treated  by  Swegle  and  Ott. 16 

Manuscript  submitted  on  September  9,  1981. 


BARKER  AND  OTTINGER 


ANOOE 


Fig.  1  —  Two  models  for  magnetically  insulated  electron  flow:  large  orbit  flow 
(left)  and  small  (armor  radius  or  Brillouin  flow  (right). 

Both  equilibrium  models  have  equal  physical  validity.  When  applied  to  a  given  problem,  they 
have  led  to  very  similar  solutions,17  however,  stability  properties  could  differ  considerably.  No 
sufficiently  extensive  analytic  comparison  has  been  performed  on  both  models  to  justify  the  exclusion 
of  one  or  the  other  for  specific  parameter  ranges.  However,  there  are  conceptual  weaknesses  regarding 
the  initiation  of  Brillouin  flow  in  coaxial  geometry.  Certainly  it  can  be  postulated  that  appropriate  injec¬ 
tion  mechanisms  can  be  experimentally  contrived,  but  this  is  not  very  satisfying  when  faced  with  the 
analysis  of  a  general,  uncontrived  device.  In  addition,  for  the  case  of  a  radial  diode  with  strong 
insulation,  the  £  x  Adrift  is  axial.  In  time,  this  would  deplete  any  initial  Brillouin  electron  population. 
Injection  and  replenishment  problems  do  not  exist  for  the  Sudan  and  Lovelace  formulation.  For  that 
reason,  it  was  chosen  for  use  in  this  analysis.  In  reality,  it  seems  probable  that  some  mixture  of  the 
two  distributions  exists  in  magnetically  insulated  coaxial  lines  with  the  Antonsen  and  Ott  picture  dom¬ 
inating  early  in  time  and  progressively  giving  way  to  Sudan  and  Lovelace  electrons  later  in  the  pulse.  A 
number  of  computer  simulations  have  been  conducted  which  seem  to  support  such  a  picture.1® 
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II.  THE  THEORETICAL  MODEL 


Cross-sections  of  both  polarity  types  of  the  coaxial  geometry  treated  in  this  analysis  are  depicted 
in  Fig.  2.  Unless  otherwise  stated,  the  central  conductor  is  assumed  to  be  the  cathode  in  the  following 
discussion.  The  analysis  for  the  case  of  a  central  anode  is  identical  except  for  magnetic  held  boundary 
conditions  discussed  at  the  end  of  this  section.  The  major  assumptions  of  the  model  are  as  follows: 


(I)  All  emitted  electrons  are  confined  to  an  azimuthally  symmetric  electron  sheath  of  radius,  rs, 
and  execute  identical  single-arc  trajectories  out  to  that  radius.  The  magnetic  insulation  field,  B, 
moderately  exceeds  the  critical  insulation  field  (defined  as  that  field  value  at  which  electron  orbits  just 
graze  the  anode  surface),  so  that  rs  <  rA.  where  rA  is  the  anode  radius  and  the  critical  fields  are 
defined  as 


and 


( B: ) 


e 


yfut  +  2  U0  , 


(®»)crii 


m0c2  0.2  a/  L/q  +  2Ua 


Here  r,  is  the  cathode  radius  and  l/„  is  the  diode  voltage  multipled  by  - — — . 

wot 


(2)  The  time-scale  is  such  that  the  self-magnetic  fields  do  not  significantly  penetrate  the  anode 
and  cathode  surfaces. 


(3)  The  electrons  are  treated  relalivistically  and  the  ions  nonrelativistically  (a  reasonable  approxi- 
matior  r'r  voltages  below  10  MV). 


(4)  The  ion  current  is  strictly  radial.  This  implies  no  significant  ion  contribution  to  the  self  mag- 

r— 

netic  field.  (The  angular  deflection  suffered  by  an  ion  is  on  the  order  of  ^  — —  radians. 

V  M()  "crii 

Therefore,  even  if  B  -  4flcm  an  ion  would  experience  less  than  0.1  radian  deflection.) 
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Fig.  2  —  Coaxial,  magnetic  insulation  geometries  of  the  negative  polarity  (top) 
and  positive  polarity  (bottom)  types. 


One  characterizes  the  electron  (ion)  beam  by  a  velocity  V  ( F),  charge  density  n(N ),  and  current 
density  j  *■=  n\(J  =  NV).  Conservation  of  relativistic  electron  energy  requires  that 


ym0c 2  -  e<f>(r)  =  m0c2. 


(1) 


where  y  «  1 1 - j 


-1/2 


.  Conservation  of  canonical  momentum  demands  that 


•y "*0v„  -  ~A„  -  0, 

where  the  II  subscript  refers  to  the  E  x  B  direction.  Finally,  conservation  of  ion  energy  states  that 


(2) 


2e 


K.  «  _  /-T7-  J<t>n  ~  <b.  (3) 

where  <f>0  is  the  imposed  diode  voltage.  Continuity  of  diode  current  flux  provides  two  additional  con¬ 
stants,  F  —  r\jr  \  and  /  -  -  r|yr|.  Employing  all  of  these  relations  allows  one  to  reduce  Poisson's  equa¬ 
tion  and  Ampere’s  law  to 


ilL 

2 

JL 

F 

r 

v. 

K 

(4) 


and 


x  x  A(r))) „  -  -  4 n 


1L 


crvr 


(5) 
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The  factors  of  two  for  the  electron  flux  terms  reflect  the  presence  of  both  radially  inward  and  radially 
outward  components  to  the  electron  current  in  the  electron  cloud.  This  doubling  of  the  electron  flow  is 
specifically  mentioned  by  Bergeron  if  Ref.  1;  however,  it  is  not  present  in  his  Eqs.  (35)  and  (36). 

Numerical  considerations  encourage  the  substitution  of  appropriately  scaled  dimensionless  vari¬ 
ables  in  place  of  those  appearing  in  Eqs.  (4)  and  (5).  Following  Bergeron,  the  terms  employed  in  this 


analysis  are  U  =  - -  d>,  W  =  - -  A„  (or  W  = - j  A .  for  an  azimuthal  insulating  field!,  and 

m0c-  m0c 2 


_  4ir  fE 
”  cB0  “V  m0 


|F|.  The  electron  flux  term  is  modified  by  the  additional  factor  of  two  mentioned 


above,  giving  K  =  — —  |/|.  A  significant  modification  to  Bergeron's  treatment  can  be  found  in  the 
cB0 

_  eB0 

new  scale  chosen  for  the  independent  radial  variable.  It  is  here  defined  as  p  =  - rr  where  B0  is 

m0cl 

chosen  to  be  some  typical  magnetic  field  strength.  In  contrast,  Bergeron  scales  r  to  the  electron  gyrora- 
dius  in  the  presence  of  the  magnetic  field  at  the  cathode  surface.  Such  a  scaling  provides  one  with  a 

dW 

simple  boundary  condition  for  the  first  derivation  of  IT  at  p(,  namely  — —  =  1.  However,  the 

value  of  B  at  p  =  p,  is  known  a  priori  only  for  the  case  of  an  azimuthal  insulation  field  with  the  cathode 
as  the  central  conductor.  In  addition,  the  possibility  of  the  magnetic  field  strength  approaching  and 
even  equaling  zero  at  the  cathode  surface  for  the  other  cases  cannot  be  discarded.  Such  an  eventuality 
would  make  nonsense  of  any  solutions  using  Bergeron's  scaling  and  may.  at  least  in  part,  explain  the 
lack  of  closed  solutions  Bergeron  noted  in  some  parameter  regimes. 


Given  these  definitions,  Eqs.  (4)  and  (5)  reduce  to 
d2U  .  dU  \(1  +  U) 


U2  +  2U  -  IT2 


U2  +  2  U 
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with  e  -  0  for  azimuthal  insulation  and  c  -  1  for  the  axial  case.  The  boundary  conditions  for  axial 
insulation  are 

-  0.  (8) 

at  the  cathode  while  at  the  anode 


U(pc)  -  W(pc) 


djU_ 

dp 


V(pA)  -  U0  and 


dU_ 

dp 


0. 


(9) 


I  Pa 


The  conditions  governing  the  first  derivative  of  U  follow  directly  from  the  assumption  of  space  charge 
limited,  Gaussian  emission  at  both  surfaces.  (However,  if  an  equilibrium  solution  is  sought  for  which 
A  =  0,  then  the  value  of  the  electric  field  at  the  anode  cannot  be  assumed  to  vanish.)  Such  an  emis¬ 
sion  mechanism  in  turn  presumes  a  sufficiently  dense  plasma  to  exist  at  both  surfaces.  We  now  have  a 
set  of  two  second  order  equations  in  U  and  W ,  with  two  additional  unknowns,  A  and  A .  Six  boundary 
conditions  are  therefore  necessary  for  a  solution.  In  addition  to  the  five  conditions  stated  in  Eqs-  (8) 
and  (9),  for  the  case  of  axial  magnetic  insulation,  the  final  condition  results  from  conservation  of  the 
magnetic  flux, 

-  Js  B  ■  da  “  §>c  <4  '  dS, 

over  a  fixed  cross  sectional  surface  area,  5,  enclosed  by  the  line  contour,  C.  Specifically,  it  is  assumed 
that  the  insulating  B-field  is  present  before  emission  begins  and  that  any  diffusion  of  the  field  into  the 
electrodes  is  insignificant  over  the  time  scale  of  the  voltage  pulse.  For  the  case  of  an  initially  uniform, 
axial  insulating  field,  B0,  which  is  taken  as  the  scaling  field  for  p,  this  condition  implies 


Pa  ~  Pc 
2  Pa 


(10) 


On  the  other  hand,  for  an  azimuthal  insulating  field,  magnetic  flux  is  not  conserved.  In  place  of 
Eq.  (10)  one  is  given  the  known  magnetic  field  strength  at  the  surface  of  the  central  conductor. 


Bcc  -  0.2  -jp. 

*V.r. 


(id 
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where  BCI.  is  in  gauss,  RC  I  is  the  radius  of  the  central  conductor  measured  in  centimeters  and  /, ,  is 
the  current  in  amperes  flowing  through  the  central  conductor.  Therefore,  the  final  boundary  condition 
for  the  azimuthal  field  case  is  simply 


dW 

df> 


B, 


(i. 

Bn 


(12) 


K'. .  °o 

where  B0  is  the  arbitrary  scale  field  and  /»,,,,  may  be  either  pA  or  pr,  depending  upon  the  geometry 
chosen.  (When  B0  -  Bt  l.  and  p,  t.  =■  pt,  this  reduces  to  Bergeron's  boundary  condition.)  For  the  case 
of  a  magnetically  self-insulated  vacuum  feed  line,  the  value  of  /„,  is  a  known  (or  at  least  independently 
derivable)  quantity.  On  the  other  hand,  /„  is  not  known  a  priori  for  a  radial  diode.  In  this  device,  the 
central  conductor  current  will  generally  be  a  function  of  the  axial  dimension,  z.  Furthermore,  the  total 
current,  ,  flowing  through  the  diode  will  be  entirely  due  to  ion  flow  through  the  A-K  gap  and  will  be 
given  in  statamps  by 


/,,  =  2n\rJ,\L. 


where  L  is  the  axial  half-length  of  the  gap.  In  order  to  clarify  the  distinctions  between  the  axial  versus 
the  azimuthal  boundary  conditions  and  unknowns,  they  are  grouped  accordingly  in  Table  1. 


For  both  cases,  the  numerical  integration  of  Eqs.  (6)  and  (7)  begins  at  the  cathode  radius,  p, 
For  a  configuration  with  B.  insulation,  values  are  guessed  for  A,  \,  and  fF'lp, ).  Then  U.  U'.  and  W 
are  integrated  out  to  the  anode  radius  and  their  boundary  values  are  compared  to  those  given  in  (9)  and 
(10).  Repeated  guesses  for  A,  \,  and  fT'lp, )  are  made  until  a  match-up  at  />A  is  achieved.  For  B„ 
insulation  with  the  cathode  as  the  central  conductor,  an  additional  starting  condition  is  given  in  the 
form  of  Eq.  (12).  This  simplifies  the  problem  considerably.  Only  the  two  anode  conditions  given  by 
Eq.  (9)  must  be  matched  and  only  values  for  A  and  A  must  be  sought  by  iterative  guesses.  Finally,  for 
the  azimuthal  case  with  the  anode  as  the  central  conductor,  the  guessed  parameters  become  A.fF'lpJ, 
and  .  The  three  conditions  which  must  be  matched  at  pA  are  now  given  by  Eqs.  (9)  and  (12).  It 
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Table  1.  Matrix  of  Knowns,  Unknowns,  and  Boundary  Conditions. 


Insulating  B-field  Type 

AXIAL 

AZIMUTHAL 

Relative  A-K  Positions 

either 

cathode  inside 

cathode  outside 

Ions  Present  ? 

Yes 

No 

Yes 

No 

Yes 

No 

Input  Parameters 

O 

f'J 

CO 

<BA 

L 

‘c.c. 

L 

Input  p=  pc  Guesses 

j,  J, 
and  W* 

j,  w 

j.  J 

j 

j,  J 
and  W 

j,  w 

Output  Quantities  to 
be  matched  at  p=p^ 

and  at  p=  ps 

U,  U\ 
and  W 

u,  u 

U,  U' 

U 

U,  U', 
and  W 

11,  w 

Universal  Input 

rft,  rc,  cp0,  LAZ,  LCIN,  and  LIONS 

Universal  Output 

Radial  Profiles  of  ip  and  B 

where:  r A  (rc)  =  anode  (cathode)  radius  in  centimeters, 

^o  =  potential  difference  in  volts, 

(Bz)o  =  imposed  insulating  axial  magnetic  field  strength  in  gauss, 

L  =  axial  length  of  A-K  gap  in  centimeters, 

Ic  c  =  current  flowing  in  the  central  conductor  in  amperes,  and 

j  (J)  =  electron  (ion)  current  density  at  cathode  in  amp/cm2. 

will  be  shown  in  the  next  section  that  analytic  solutions  are  available  for  the  region  between  and  ps. 
This  allows  for  numerical  solution  matching  at  ps  instead  of  at  i>A  with  a  considerable  gain  in  accuracy. 

A  problem  arises  when  solutions  are  sought  for  radial  diodes  in  which  both  axial  and  azimuthal 
magnetic  insulating  fields  are  present.  If  no  B.  is  imposed  in  the  A-K  gap.  then  only  pure  azimuthal 
insulation  is  possible  and  no  conflict  can  arise.  When  an  axial  insulation  case  with  ion  flow  is  being  stu¬ 
died,  however,  it  is  inevitable  that  a  B„  be  generated  as  well.  The  strength  of  this  BH  in  the  A-K  gap 
can  be  predicted  directly  from  the  value  of  \  obtained  in  the  pure  B.  solution.  If  this  field  strength  is 
comparable  to  or  greater  than  the  imposed  B:.  the  numerical  formulation  presented  here  has  more  lim¬ 
ited  value  for  predicting  actual  diode  operating  characteristics.  (See  Appendix  B.) 
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III.  THE  COMPUTATIONAL  FORMULATION 

As  previously  stated,  the  solution  of  Eqs.  (6)  and  (7),  combined  with  the  boundary  conditions  in 
Eqs.  (8)  and  (9),  will  be  carried  out  as  a  three-parameter  "shooting"  problem.  Values  must  be  guessed 
dW 

lor  A,  A,  and  and  the  equations  then  integrated  from  p,  to  pA.  The  correctness  of  a  given 

guess  is  determined  by  the  accuracy  with  which  generated  function  values  match  with  known  solutions 
at  pA  and  ps.  In  undertaking  the  numerical  integration  of  Eqs.  (6)  and  (7),  special  attention  must  be 
paid  to  three  separate  values  of  p  at  which  certain  terms  on  the  right  hand  side  become  singular. 
Specifically,  the  denominator  of  the  electron  current  term  vanishes  at  p  =  p, ,  since  both  U  and  W  are 
defined  to  be  equal  to  zero  there.  It  also  vanishes  at  the  electron  sheath  boundary,  p  =  pv,  as  would  be 
expected,  since  this  quantity  represents  the  radial  component  of  the  electron  velocity.  The  final  singu¬ 
larity  occurs  at  p  *=  pA  where  the  ion  current  it/  t  in  Eq.  (6)  becomes  infinite.  Brute  force  numerical 


procedures  cannot  cope  with  integrations  in  the  neighborhood  of  these  singular  points.  Instead,  analytic 
treatments  expanding  about  the  singularities  produce  approximate  explicit  solutions  valid  in  the  patho¬ 
logical  regions  nearp,  andp^.  This  technique  provides  the  numerical  quadrature  with  tractable  starting 
values  at  p  =  p,  +  Ap.  The  numerics  may  be  relied  upon  from  that  point  up  to  about  p  =  pv  -  Ap  at 
which  radius  another  analytic  approximation  advances  the  solution  to  px.  Beyond  that  point.  A  equals 
zero.  Equations  16)  and  (7)  have  straightforward  analytic  solutions  spanning  pv  <  p  <  pA  which  can 

be  used  to  accurately  advance  U ,  W  or  — —  ,  and  at  p  =  p<  out  to  the  anode  radius.  In  this 

l  dp  dp 

manner,  the  numerical  singularity  at  p  =  pA  need  never  be  faced  directly. 


The  numerical  solution  of  Eqs.  (6)  and  (7)  begins  with  the  derivation  of  an  explicit  functional 


expression  for  the  radial  variations  of  U  and  W  in  the  immediate  vicinity  of  the  cathode  surface.  For 
p  “  p,  +  r  where  r  «  p(  ,  the  equations  can  be  expanded  and  solved  in  terms  of  a  triple  power  series 


with  solutions 


'no  ,  n  oo  »  s/3  „  n 

Id-  f-  X"."1  -  +  -  I  -  +  —  I  U„,5/3'  . 

P‘  n- 0  Pc  Pc  I  (i-0  Pc  Pc  „_o  Pc 
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„  n  14/3  M  (  n  5/3  „ 

W - -  Ti,"’  -  +  —  T  i,w,)  -  +  —  Tt>J5/}>  — 

Pc  n-0  Pc  Pc  n-o  Pr  Pc  n-0  Pc 


Straightforward  calculus  and  algebra  may  then  be  used  to  isolate  the  leading  order  terms. 


i/-  4 
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where  the  term  bc  =  is  presently  unknown  and  where  the  solutions  are  only  valid  for 
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9A 
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3 
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Equation  (18)  places  an  upper  limit  on  the  initial  step  size.  All  but  the  ion  term  in  Eq.  (16)  had 
been  predicted  by  the  leading  order  terms  of  an  expression  due  to  Goldstein19  which  corrected  the 
Child-Langmuir  emission  law  to  include  transverse  magnetic  field  effects.  The  first  derivatives  of  equa¬ 
tions  (16)  and  (17)  are  simply 


1  JL 

1/3 

aj/3  + 

r 

4  Pc 

10 

PcU^\ 


w.  _  dW  dW 
dp  dr 


“(-D**1*, i--  +  T  hw  - 

dr  p,  4  9  p, 


For  purposes  of  the  overall  numerical  quadrature  of  Eqs.  (6)  and  (7),  the  gap  region  between  p,  and 
pA  is  divided  up  into  some  predetermined  number  of  identical,  one-dimensional  cells  of  width,  Ap. 
The  above  approximations  are  used  to  advance  the  solution  of  U,  W,  U\  and  W'  over  the  first  Ap  away 
from  the  p  =  pc  singularity  point. 

The  bulk  of  the  remaining  numerics  is  accomplished  using  a  proven  predictor-corrector  scheme. 
Unfortunately,  this  scheme  requires  a  knowledge  of  functional  values  at  four  prior  spatial  steps  in  order 
to  advance  an  additional  Ap.  In  order  to  obtain  those  additional  values  at  the  start  of  the  integration,  a 
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fourth-order  Runga  Kutta  scheme  is  employed.20  For  this,  as  well  as  for  the  predictor-corrector,  the  set 
of  two  second-order  D.E.’s  is  reduced  to  a  set  of  four  first-order  D.E.'s  as  follows: 


U 1  s  X, 


v,  K(  1  +  U)  \  X 

pVs  p  • 

W'  =  Y.  and  (23) 

+  (24, 
f)\/S  f}J  ft 

where  S  =  U2  +2  U  —  W'2.  Each  expression  is  then  in  the  form  ry’  —  /,(/*, U.  W.X,  Y)  and  can  be 
advanced  one  Ap  at  a  time  using  the  standard  R-K  approximation. 


n»+i  -  +  y<*,,i  +  2At)2  +  2A:^3  +  A,*). 


where 


A„,  =  lp/v{pH,U„.tV„,Xn,Y„), 

A„2  =  V/„(p„  +  +  y*ti.  W„  +  y-A„|.  ...), 

=  W„<P„  +  y Xp.Vn  +  yAf2.  w„  +  yArir.v  •  •). 

*r,4  ==  Ap./^p,,  +  Ap.f/„  +  kV),  Wn  +  Avj,  . . .). 

Thereafter,  the  Adams-Bashforth-Moulton  predictor-corrector  scheme21  advances  the  variables  in  a 
iwo-step  process. 

(^/i+lVedicied  =  V„  +  y^(5-<'r?„  ~  59t}„_,  +  37t)„.2  -  9r),' 


An 

^n  +  |)correcied  =  »»«  +  yy  (9'9«  +  l  +  ~  5rj' 

where  tj„+i  has  been  calculated  using  the  predicted  value  for  i7„+|.  An  error  check  is  accomplished  by 
monitoring  the  magnitude  of  \(y)p  -  n,  )lr)n\  where  r)p  and  77,  are  the  predicted  and  corrected  values  of 
Tj,  respectively. 


The  above  predictor-corrector  process  is  continued  until  the  electron  sheath  radius,  pv,  is  reached. 
At  that  point,  another  analytic  approximation  must  be  found  in  order  to  correctly  advance  the  function 
values  past  the  numerical  singularity  there.  This  is  accomplished  by  expanding  Eqs.  (6)  and  (7)  about 
(>s  'r-  a  power  series.  One  may  write  p  =  ps  -  r,  U  «  £/s  +  0,  and  W  -  +  W  where 


II 
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Ps  »  r  >  0, 

0  -  a0ra°  +  0|f“‘  +  a2r“2  +  . . .  , 
W  -  6</°  +  b/'  +  +  ...  , 


and 


Uj  +  2f/s  -  =  0. 

Equations  (6)  and  (7)  can  then  be  expanded  and  solved.  The  leading  terms  of  the  solutions  are 


U  **  Us  +  a0r  +  axr 3,2  +  air1. 


and 


where 


and 


W-  Ws  +  bQr  +  bxrm  +  b2r2. 

2(1  +  Us)a0-  2Wsb0  >  0. 

4x(l  +  Us) 
a  i 

bi 

a2  = 
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dS) 
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(a0  —  \/(Up—  Us)'n) 
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( b0  -  c  ITs/ps) 
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As  is  obvious  from  Eqs.  (25)  and  (26),  the  coefficients  a0  and  b0  are  merely  the  negative  first  deriva¬ 
tives  of  U  and  W,  respectively  at  ps  (  i.e.,  a0  =  -  £/$  *=  -  Xs  and  b0  =  -  —  -  Ks).  Having  used 

Eqs.  (25)  and  (26)  to  obtain  and  Ks,  the  numerical  integration  of  Eqs.  (6)  and  (7),  with  K 

now  set  equal  to  zero,  may  resume.  In  the  RADBER  code,  however,  explicit  analytic  expressions  are 

used  to  match  all  functions  of  interest  at  ps  from  the  known  values  at  the  anode  radius. 


Beyond  the  electron  sheath  radius,  the  key  deterministic  equations  for  the  electrostatic  and  mag¬ 


netostatic  potentials  become 
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Equation  (27)  takes  the  same  form  as  the  equation  treated  by  Langmuir  in  his  1913  paper23  and  later 
treated  in  greater  depth  for  coaxial  geometry  in  a  subsequent  paper  by  Langmuir  and  Blodgett.23  The 
solution  to  the  equation  is  therefore 

( 9  \in 

l/=(/0-  jW 32A  (29) 

where  j3  =  L/Ly"  and  y  =  In  MM.  The  first  14  values  of  the  coefficients,  An ,  are  listed  in  Blodgett's 

M 

paper.  Knowing  the  values  of  t/„,  A,  pA ,  and  ps,  a  value  for  Us  may  be  calculated  using  Eq.  (29). 
This  value  is  then  compared  to  that  found  via  the  numerical  integration  of  Eq.  (6)  out  from  p, .  The 
closeness  of  the  match  is  used  to  measure  the  accuracy  of  the  original  set  of  parameter  guesses.  For 
the  case  where  \  =  0  (i.e.,  no  ion  flow)  Eq.  (27)  becomes 


with  solution 


P 


(30) 


u  -  U5  +  xsps  In  M- 
J  P.v 

and  the  same  matching  procedure  for  Us  may  be  carried  out. 
the  form 


(31) 

Similarly,  the  solutions  to  Eq.  (28)  take 


W  = 


M's- 


<?[)! 

PS 


JL  +  £»  e  =  I 


Ws  +  tf„2  In 


+ 

P.v  P 


Ps 


(32) 


e  =  0 


Matching  the  solutions  for  the  first  derivatives  at  p  =  ps  demands  that  a oi  =  y  (  ^tP.v  ~  Fsp.v)  for  the 


*  =  1  case  and  a02  =  TsP.v  (°r  the  «  =  0  case.  Equations  (29),  (31),  and  (32)  may  be  used  to  obtain 
values  for  Us,  Xs.  Ws,  and  Kv  in  terms  of  known  quantities.  Depending  upon  the  specific 
configuration  under  study,  some  subset  of  those  four  values  (see  Table  I)  are  matched  to  the 
corresponding  values  on  the  cathode  side  of  the  electron  sheath  boundary  as  obtained  from  Eqs.  (25) 
and  (26).  The  closeness  of  that  fit  at  p  -  ps  of  the  two  sets  of  solutions  forms  the  final,  most  accurate 
lest  of  the  validity  of  the  original  parameter  guesses. 
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IV.  APPROXIMATE  ANALYTIC  SOLUTION 


For  the  case  of  low  voltages  and  with  geometries  whose  anode-cathode  gap  is  much  smaller  than 
the  radius  of  the  outer  conductor,  complete,  approximate  analytic  solutions  can  be  obtained  by  match¬ 
ing  the  analytic  solutions  near  the  cathode  given  in  Eqs.  (16M20)  with  the  solutions  valid  near  the 
electron  sheath  radius  given  in  Eqs.  (25)  i...i  (26).  The  solutions  are  matched  at  pt  —  (p,  +  ps)/ 2  by 
requiring  that  U ,  W,  U\  and  B/B0  =  ( W'  +  «  W/p)  all  be  continuous  at  p  -  p|.  At  p  =  Ps  the  solu¬ 
tion  given  in  Eqs.  (25)  and  (26)  is  also  matched  in  the  same  manner  with  the  solution  valid  for  ps  ^ 
p  <  pA.  There  one  finds 


I9  * 

T 

2/3 

4/3 

U-U0- 

r 

(33) 

1  ^ 

Pa 

and 


W  = 


M's  +  p  *  1° 


for  e  =*  0 


(34) 


where  r  =  pA  -  p  and  bA  =  B(pA)/B0.  Equation  (33)  is  equivalent  to  the  limit  of  Eq.  (29)  for  r  « 
pA.  Likewise,  Eq.  (34)  follows  from  Eq.  (32)  when  a0  is  obtained  by  imposing  the  boundary  conditions 
at  p  =  pA  instead  of  at  p  -  ps.  For  «  =  1,  flux  conservation  fixes  the  value  of  )V(p  =  pA)  (see  Eq. 
(10))  and  yields  bA  -  [(pj  -  p})/(pA  -  pj)  -  2psWs/(pl  -  pj)].  For  the  azimuthal  field  case  («  = 
0),  on  the  other  hand,  bA  =  YA  in  agreement  with  Eq.  (32)  since  the  field  outside  the  sheath  is  con¬ 
stant.  In  order  to  complete  the  overall  solution  for  this  analytic  approximation  one  must  also  use  the 
relationship,  L/£  +  2  Us  -  H's  =  0.  and  also  bc  =  1  for  the  e  -  0  case.  There  is  now  enough  informa¬ 
tion  with  which  to  construct  a  composite  solution. 


For  the  case  where  the  central  conductor  is  the  cathode,  one  finds  that  for  both  t  =  0  or  1 ,  an 
analytic  solution  is  possible  as  long  as  the  diamagnetic  electron  current  does  not  significantly  alter  the 
magnetic  field.  In  this  case  B0  is  approximately  given  by  the  applied  field. 
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J(pa') 
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|J/2 
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where  d  =  pA  -  pc  and  jCL  =  — 


m, 


ps-pc  =  (Us2  +  2Us)'-2. 
1/2 


(36) 

(37) 

(38) 


m. 


JCL  =  (2e/me)ui  d>ll2l<in(rA  -  r,  )2.  The  sheath  potential 


needed  to  complete  Eqs.  (36M38)  is  defined  by  the  transcendental  equation 

6  d 


Us 

U°=-f 


I  + 


(39) 


(Us2  +  2Us)'n 

The  ion  production  efficiency  of  the  diode,  rj  =  J/U  +  j),  can  be  easily  obtained  Trom  Eqs.  (36)  and 
(37).  Plots  of  J  and  ps  —  p(.  are  shown  in  Figs.  (3)  and  (4)  as  solid  lines.  The  validity  of  the  results 
becomes  suspect  beyond  e<t>o/ m,.c2  ~  1.  Note  that  Ip A  —  p,  |  =  d  «  pA  was  also  assumed  which 
implies  these  results  only  apply  to  large  radius  structures  with  small  anode-cathode  gaps  (i.e.,  large 
aspect  ratio  devices).  For  higher  voltages  and  smaller  aspect  ratio  devices  the  RADBER  code  is  needed 
since  in  these  cases  the  simple  matching  technique  used  to  form  the  composite  analytic  solution  is  no 
longer  possible.  This  results  from  the  decreasing  range  of  validity  of  the  expansions  used  to  obtain  the 
solution.  In  any  case,  however,  these  approximate  solutions  may  be  used  to  provide  reasonable  starting 
parameter  guesses  for  the  RADBER  code. 
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t'ig.  3  —  Analytic  estimates  for  ion  current  enhancement 
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Fig.  4  —  Analytic  estimates  for  electron  sheath  thickness 


V.  THE  NUMERICAL  TREATMENT 

The  RADBER  computational  package  consists  of  a  main  core  program  and  four  subprograms: 
START,  PRIME,  BOUND,  and  MATCH.  All  of  these  program  elements  are  listed  in  their  entirety  in 
Appendix  A.  They  are  written  in  a  version  of  FORTRAN  compatible  with  the  Texas  Instruments 
Advanced  Scientific  Computer.24  Every  effort  was  made  to  avoid  non-standard  terminology.  Neverthe¬ 
less,  the  need  for  minor  changes  can  be  anticipated  when  implementing  this  code  at  any  other  computer 
facility.  Particularly  prone  to  such  changes  are  the  data  input/output  statements.  For  that  reason,  I/O 
formats  were  left  as  simple  as  possible.  A  particularly  strong  point  of  the  code  is  its  vectorized  format. 
This  allows  for  the  simultaneous,  step-by-step  integration  of  U  and  W  for  a  large  number  of  X  and  A 
guesses  in  a  "shotgun"  fashion.  The  average  cost  for  the  complete  integration  of  Eqs.  (6)  and  (7)  over 
500  steps  from  pc  to  pA  for  a  single  (X,  A,  W'(pc))  guess  was  about  $0,025.  (CPU  rates  are  over 
$900/hour  on  the  ASC  at  NRL.) 

The  basic  flow  chart  for  RADBER  is  shown  in  Fig.  5.  In  the  main  program,  data  input  begins 
with  a  choice  for  the  parameter,  NDR,  the  number  of  radial  steps  across  the  anode-cathode  gap.  The 
default  value,  as  shown,  is  500.  The  specific  physical  and  numerical  parameters  for  the  problem  under 
study  are  fed  in  through  the  NAMELIST’s,  GUESS  and  GIVEN.  For  the  sake  of  practicality,  the  code 
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fig  5  —  RADHER  flow  chart 


accepts  guesses  for  the  actual  current  densities,  ./  and  J ,  at  the  cathode  rather  than  for  the  dimension¬ 
less  fluxes,  X  and  A.  In  GUESS,  the  range  of  guesses  for  the  triplets,  (./,  7,  H-"(p,))  are  specified.  The 
minimum  guess  for  the  electron  current  density  is  AMIN  in  amperes  per  cm.2  A  total  of  NA  values  of  ./ 
are  then  tried  in  the  program,  each  one  DELA  greater  than  the  next.  Similarly,  there  are  NB  guesses 
for  J  in  amps/cm2  beginning  with  BMIN  and  increasing  in  steps  of  DELB.  The  product  NA  multiplied 
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by  NB  must  equal  exactly  100  for  the  program  as  presently  written.  Finally,  the  slopes  of  W  at  the 
cathode  radius  are  guessed  as  some  fraction,  FMIN,  FMIN  +  DELF,  FMIN  +  2  DELF,  etc.,  of  the 
scale  magnetic  field,  B0.  There  is  no  limit  on  the  number,  NF,  of  these  which  can  be  tested  with  the 
100  (J.J)  guesses  under  consideration.  The  significance  of  the  integers,  11,  12,  and  13  will  be  explained 
when  the  OUTPUT  section  of  the  program  is  discussed. 

In  NAMELIST,  GIVEN,  the  cathode  and  anode  radii  are  specified  in  centimeters  by  RC  and 
RANODE,  respectively.  The  potential  difference  between  the  two  electrodes  is  given  in  volts  by  U0. 
The  scale  magnetic  field  strength,  Ba,  in  gauss,  is  arbitrary,  but  can  conveniently  be  chosen  as  the 
vacuum,  imposed  field  strength  for  an  axial  insulation  field  problem  or  as  the  vacuum  value  near  the 
cathode  for  an  azimuthal  insulating  field.  The  parameter,  GAPLEN,  is  only  significant  when  azimuthal 
insulation  in  a  radial  diode  configuration  (i.e.,  the  cathode  is  the  central  conductor)  is  under  study.  In 
that  case,  it  is  equal  to  the  axial  length  of  the  A-K  gap  in  centimeters.  In  a  similar  sense,  CURCEN 
need  only  be  specified  for  the  case  of  azimuthal  insulation  with  no  ion  flow  present.  It  is  equal  to  the 
axial  current  assumed  to  be  present  in  the  central  conductor  in  amperes.  The  logical  variables,  LAZ, 
LCIN,  and  LIONS  specify  the  overall  nature  of  the  problem  under  consideration  as  follows: 

a.  LAZ  —  TRUE,  implies  an  azimuthal  magnetic  insulation  field, 

b.  LCIN  -  .TRUE,  implies  the  cathode  is  the  central  conductor,  and 

c.  LIONS  *  .TRUE,  implies  the  presence  of  ion  (proton)  flow. 

Finally,  ERROR  is  a  decimal  fraction  to  be  chosen  as  the  maximum  allowable  percentage  difference 
tolerable  in  predictor-corrector  steps  before  a  given  step  iteration  there  is  stopped. 

RADBER  then  uses  the  above  NAMELIST  inputs  in  order  to  calculate  constants  and  parameters 
utilized  throughout  the  rest  of  the  program.  These  set  of  one  hundred  guesses  for  (J.J)  are  then 

dW 

attempted  for  each  guess  of  — — (p  -  pr).  That  is,  the  entire  remainder  of  the  program  is  cycled 

dp 

through  once  for  each  value  of  W'Tpf).  This  iteration  begins  at  statement  number  400.  The  100 
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integrations  from  pc  to  pA  for  each  ( j,J )  pair  are  carried  out  simultaneously.  Quadrature  over  the  first 
Ap  is  accomplished  in  the  subroutine,  START,  using  the  analytic  approximations  of  Eqs.  (I6M20). 
From  that  point  and  over  the  next  4  Ap,  fourth-order  Runga-Kutta  is  employed  to  generate  a  finite 
length  string  of  starting  values  for  the  master  predictor-corrector  scheme.  Both  here,  as  well  as  in  the 
main  integrator,  the  values  for  the  second  derivatives  of  U(p)  and  W(f>)  are  calculated  via  Eqs.  (6) 
and  (7)  in  the  subroutine,  PRIME.  That  subroutine  also  locates  the  electron  sheath  radius,  ps,  by 
monitoring  the  value  of  U2  +  2U  -  IF2  and  finds  the  effective  anode  radius  for  a  given  lj,J)  by  pin¬ 
pointing  the  value  of  p  for  which  U(p)  =  Un.  The  Adams-Bashforth-Moulton  predictor-corrector  algo¬ 
rithm  integrates  "over"  p  =  ps,  ignoring  any  functional  discontinuities  there.  It  simply  forges  ahead  in 
a  brute  force  manner,  setting  A  equal  to  zero  between  pv  and  />A.  For  some  choices  of  parameter 
ranges,  this  may  lead  to  serious  errors  propagating  in  the  quadrature  beyond  p  =  py.  For  that  reason, 
the  array,  IFLAG,  stores  the  integer  number  of  the  radial  step  over  which  the  electron  sheath  boundary 
was  encountered.  The  final  test  for  the  validity  of  the  various  (/,/  H/'(pi ))  guesses  is  accomplished 
through  matching  of  the  potentials  and  the  fields  at  p  =  ps. 
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After  the  predictor-corrector  has  blindly  pushed  all  100  integrations  out  to  p  =  p4,  the  program 
looks  back  to  find  the  function  values  at  the  two  data  points  immediately  preceeding  each  of  the  100 
respective  py’s.  Having  stored  those  values  in  the  COMMON  block,  FIT,  the  program  then  calls  the 

P 

subroutine,  BOUND,  to  calculate  V,  U'(=X),  W,  and  H /,(=Y)  at  p.s  (approaching  from  the  cathode  j| 

side)  via  the  analytic  approximations  given  in  Eqs.  (25)  and  (26).  Specifically,  the  lowest  order  terms  - 

of  those  equations  may  be  rearranged  to  yield  ■ 

I 

Us  =  +  Xsr  +  a,  r'/2,  (40) 

Ws  =  Wn +  Ysr  +  rV2.  (41)  i; 

t- 

where  the  subscript  (/;  -  I),  signifies  the  function  value  at  the  second  data  point  before  the  recorded  j. 

py.  Equation  (40)  is  multiplied  by  By  and  (41)  by  (I  +  Us ) ■  Subtraction  of  the  resultant  equations  ,< 

leads  to  the  expression 

^-,(1  +  Us)  ~  j»s<l  +  (42)  ' 

Ty  (Fy  -  Fy(l  +  Uy) 
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but,  since  U£  +  2US  -  -  0, 

,  _  Wi  +  2Us)''i{\  +  Um_t)-  K-id  +  UJ 

vsu  +  us)  -  xs(us2  +  2vs)l/2  3 

Using  the  known  values  of  U„  as  reasonable  first  guesses  for  Us,  Eq.  (43)  yields  a  first  predicted  value 
for  r.  That  value  can  then  be  fed  into  (25)  to  obtain  a  corrected  value  for  Us  and  the  process  contin¬ 
ued  until  reasonable  convergence  is  achieved.  As  written,  the  subroutine  BOUND  performs  a  number 
of  iterations  equal  to  the  integer,  ITER,  which  is  a  fixed  parameter  in  the  CALL  statement. 

At  that  point  in  the  RADBER  program,  a  pause  is  made  in  the  computations  in  order  to  output 
the  results  of  the  full  integration  to  p  =•  pA.  The  format  of  this  primary  output  depends  upon  the 
parameters  and  geometry  under  consideration  as  indicated  in  Table  1.  In  all  cases  the  guessed  values 
for  j  and  J  (A  and  B)  in  amps/cm2  are  listed  as  are  the  values  found  for  rs  and  rA  in  centimeters  as 
well  as  the  values  of  the  electric  potential  and  electric  field  at  p  =  pA  in  volts  and  volts/cm  respectively. 
Depending  upon  the  geometry  being  studied,  W(pA)  or  tV'(pA)  may  be  listed  along  with  the  value  it 
must  have  in  a  true  steady  state.  The  table  thus  created  allows  a  quick  evaluation  of  the  merits  of  the 
(J,J)  guesses  tried  for  the  given  value  of  W"(pt).  In  addition,  immediately  following  that  table  are 
listed  complete  radial  profiles  for  the  electric  potential  in  volts  and  for  the  insulating  magnetic  field  in 
gauss  from  pc  to  pA  for  the  1 1th,  12th,  and  I3th  < j.J )  guesses.  These  profiles  provide  invaluable  infor¬ 
mation  about  how  well-behaved  the  steady  slate  conditions  are  for  typical  current  density  choices. 

Having  completed  those  pieces  of  data  output,  RADBER  moves  on  to  a  comparison  of  the  func¬ 
tion  values  at  the  sheath  boundary  obtained  by  numerical  integration  on  the  cathode  side  and  by  ana¬ 
lytic  evaluation  in  the  electron-free  region  on  the  anode  side.  The  analytically  obtained  values  are  cal¬ 
culated  in  the  subroutine,  MATCH,  using  Eqs.  (29)-(32)  along  with  the  values  for  ps,  Xs,  and  Ys 
obtained  in  the  subroutine,  BOUND.  These  pairs  of  values  to  be  matched  at  the  sheath  boundary  are 
tabulated  as  the  final  block  of  output  in  a  given  H/,(p, )  cycle.  The  table  heading  clearly  indicates  which 
values  should  be  compared  to  which  (e.g.,  US l  to  US2).  It  is  useful  to  use  this  sheath  boundary 
matching  to  fine-tune  (j,J,W’(pc))  guesses  after  some  coarse  fit  to  the  chosen  parameters  has  been 
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established  in  the  p  =  i>A  matching  table.  The  process  suggested  for  this  iterative  improvement  of 
guesses  is  illustrated  by  the  examples  given  in  Appendix  B  and  Appendix  C. 

VI.  CONCLUSIONS 

The  RADBER  computer  code  which  has  resulted  from  the  extensive  theoretical,  computational, 
and  numerical  efforts  outlined  above  stands  as  a  unique  tool  for  analyzing  numerous  practical,  high- 
voltage  devices.  It  provides  the  experimenter  with  reliable  predictions  for  electron  and  ion  flow  profiles 
as  well  as  electric  and  magnetic  field  characteristics  in  coaxial,  cylindrical  geometry.  Appendix  A  which 
follows  this  section  provides  a  complete  FORTRAN  listing  of  the  code.  Appendix  B  illustrates  the 
details  of  how  the  code  is  applied  to  a  specific  problem  and  how  the  iterative  "shooting"  technique  is 
used  to  converge  to  the  correct  answer.  The  device  chosen  in  this  first  sample  case  is  a  simple  radial 
diode.  Its  treatment  points  out  the  danger  of  applying  RADBER  to  a  mixed  axial-azimuthal  insulation 
field.  Finally,  the  much  simpler  problem  of  electron  flow  in  a  vacuum  transmission  line  is  treated  for 
three  different  voltages  in  Appendix  C.  The  solution  has  direct  application  to  a  power  loss  problem 
experienced  in  a  more  complicated  coaxial  geometry.  Taken  together,  the  two  sample  treatments 
explicitly  demonstrate  the  versatility  of  the  code. 
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THE  RADBER  COMPUTER  CODE 

C 

c 

- - 

C  THE  RADBER  PROGRAM 

C 

PARAMETER  NDR-750 
C 

LOGICAL  LAZ.LC1N, LIONS 
INTEGER  IFLAGdOO) 

REALXS  R, DR, RC, RANODE 
C 

DIMENSION  X(  100, NDR >,V(  1O0.NDR ),XP(  100.NDR >,  VR<  100, NOR ),URA<1O0> 
k  .uueo.NCRi.uuoo.NDRi.uudooi.uudOOJ.xxdooi.YVdeo) 

I  ,AC(103>,BC<100>,  SCAd00),SC8d00  .SCCd0e>,SCDd00>, 
k  SCEdCO).SCF<100),SCGd00>,SCHd0O> 

c 

COMMON  /CCEF  /  AC100  ) ,  B<100  ),RS(I0O  >,RA(10O >,  AFLACdOO), 
k  R , DR , RC , RANODE , R I MU . UO . B0 , EBVMC2 . EHRAT.C . GAPLEN 

COMMON  /SUITCH/  LA2, LCIN, LIONS 

COMMON  /FIT/  UN(  103  )  ,UNdC3  ),  XNdOO  > » VM(  160 >,UNfll ( 100  ),UNH1 ( 100), 

&  XHM 1 ( 1 00 ) , VNM1 1 1 00 ) , RNM1 d  00 ) , RST AR  d  00 ) , USTAR (100), USTAR ( 100 ), 

k  XSTARt 100),VSTOR< 100 ) 

COMMON  /ITCHY/'  SCRA<  100  ).US(10O  ),US<100  ),ES(100  ) 

c 

NAMELIST  /GUESS/  AMIN. DELA, NA, BMIN, DELB,NB.FMIN,DELF,NF. II , 18, 13 
NAMELIST/GIUEN/RC.RANODE.U0, B0, GAPLEN, CURCEN, ERROR, LAZ, LCIN, LIONS 
C 

C  IN  NAMELIST  'GIUEN'  - 

C  GAPLEN.  RC.  AND  RANODE  ARE  IN  CENTIMETERS, 

C  UO  IS  IN  COLTS, 

C  BO  IS  IN  GAUSS. 

C  CUPCEN  IS  IN  AMPERES, 

C  ERRCR  IS  A  DECIMAL  FRACTION.  AND 

C  LAZ-.TRUE.  FOR  AN  AZIMUTHAL  INSULATION  FIELD. 

C  LCIH* .TRUE .  FOR  THE  CATHODE  AS  THE  CENTRAL  CONDUCTOR 

C  LIONS*. TRUE.  FOR  THE  FRESENCE  OF  ION  FLOU 

C 

c 

C* ttttttttt  SUPPLY  NEEDED  CONSTANTS  ************************* *1 

C 

sixinv- 1.00/e. 00 

C-2.9979E10 
PI*3.  1416E0 

EMRAT- SORT ( 1 . S72GE-84/ <8.0*9. 1095E-88  )  ) 
£BVMC2-4.8032E-10/<g.l095E-88*C*C> 

C 

Ctttt*****  INITIALIZE  AND  SCALE  THE  VARIABLES  XtXXXXttXXXXttTXXXXXXXX 
C 

CXXX  FIRST  READ  IN  THE  DEFAULT  PARAMETERS  *** 

C 

READIS. GUESS) 

'  REAOCS, GIVEN) 

C 

C*M  THEN  SUPPLY  THE  SPECIFIC  VALUES  *** 

C 

READ (S, GUESS) 

READ<S, GIVEN) 

ASCALE  •  8.0  *  PI  *  RC  *  3.0E9  /  (C  *  BO) 

BSC ALE  •  ASCALE  *  0.S  *  EMRAT 
RANODE  •  EBVr.C<2*30*RAN0DE 
RC  •  E8YMC2X8CXSC 
U3  •  EBv;X£*UO/3.0E2 
DR-lRF.NOIE-RC  )/DFLOAT(NDR) 

FACT-DR/24.0 
URITEC6, GUESS) 

URITECS, GIVEN) 

C 

tmsCENERATE  THE  CURRENT  FLUX  COEFFICIENTS  ********************* 

C 

DO  2  I-l.HA 

AA-  <AMIN+DELA*< 1-1 )  )*ASCALE 
DO  8  J-l.NB 
AC(I+NB*(J-1))  •  AA 

BC(I*NB*(J-1))  -(BMIN+DELB*<J-1)  )*BSCAL£ 

8  CONTINUE 
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CXXXXXXXXXX  CALCULATE  CORPECT  U<RA>  OR  IIPIRA)  lUttttUWUW* 

C 

IF ( LAZ )  CO  TO  500 

Ur  AC- (RAN0DEXX2-RCXX2  )/(2.0XRANODE ) 

CO  TO  50 t 

500  IF(LCIN)  GO  TO  503 
IF<LI0NS)  GO  TO  503 

UF AC-3040. 2XCURCENXEBYHC2/RAN0DE 

501  DO  504  1-1,100 

504  U3A( 1  )-L'FAC 
GO  TO  se2 

503  UF AC -0 . 0SXEBYMC2XCXB0XCAPLEN/ ( EHRATXRANODE ) 

DO  SOS  1-1,100 

505  U?A< I )-UFACXB<  1 ) 

502  CONTINUE 
C 

cmmtm  iterate  ouer  choices  for  cathode  b-field  iimmmmm 

c 

IT-0 

400  IT-IT-H 

F-FHIN*D£LFX<IT-1) 

R-RC 

IF < LAZ. AND . LC IN . AND . . NOT . LIONS )  F-0 . 2XCURCENXEBYHC2/RC 

ICORR-0 

DO  1  1-1,100 

RS(I  )*0.0 

RACD-0.0 

AFLAGa>-1.0 

IFLAGU  )-0 

A( I )-AC( I ) 

B(I)-BCd) 

1  CONTINUE 
C 

CXXXXXXXXXX  U5E  ANALYTIC  APPROX I HAT IONS  TO  CENERATE  INITIAL  XXXXXXXXXXXX 

cxxxxxxxxx  ualues  at  rc+dr  mmttmmmmmmmttt 

c 

R-R+DR 

RINU-1.0/R 

C 

CALL  START<U,X,U,Y,F) 

C 

CXXXXXXXXX  TH  ORDER  RUNGAHCUTTA  FOR  STARTING  UALUES  XXXXXXXXXXXXXXXX 
C 

DO  4  1-1,4 
C 

DO  5  J-1.I00 
SCA(J)-DR*X(J,I) 

5  SC8( J )-DRXY( J, I ) 

DO  6  J- 1,100 

SCC< J  )-DRX(X( J, I )+SCA( J 1X0,5 ) 

6  SCD<  J )-DPX( Y( J, I )+SC9( J  )X0.S ) 

DO  7  J-1,103 

SCE(J)-DR*<X(J,I)40.5XSCC(J) J 

7  SCF ( J  )-DRX(Y( J,I 1+0.5XSCDCJ) > 

DO  8  J* I , 100 

SCG<  J  l*DRX(X(  J,I HSCE( J )  1 

8  SCH<J)-DRX<V(J,I)4SCF(J>) 

c 

DO  9  J-1,100 

U1J.I+1  )*U(J,I)»(SCA(J)+a.0XSCC<J)*a.0XSCE(J)4SCG(J>)XSIXINU 

9  U(J,I+l)-U(J,I  )*<5CB( J  )+2.0XSCD(  J  >+2.0XSCF(  JHSCHC  J>  1XSIXINU 
G 

CALL  PRinE(U<l,I),U(l,I)>X(l,n(V(l,I),XP(l#I)#YP(l,I)»t#IFLAG) 

c 

DO  10  J-1,100 
SCGC J )*DRXXR< J, I ) 

10  SCH( J  )-DRXYP< J, I ) 

C 

R  -R+0.5XDR 
RINU-1.0/R 
C 

DO  11  J-1,100 
UU(J)-U(J,I )+0.SXSCA( J ) 
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C 

C 


16 

C 

4 

C 


UU<J>-UU,I)+O.StSCI(J> 

XX  ( J  >-X  JL  : )40.SISCG( J ) 

vv(j)-YtJ,i)*e.sxscH(j) 

CALL  PRXNE<UU,UU,XX,YY,SCA,$CI,I,IFLAG) 

DO  18  J-l, 10* 

SCA( J  )-DRXSCA( J ) 

SC«(J)-DRXSC0(J> 

DO  13  J-l, 100 

UU( J )-U< J, I )49.5XSCC( J) 

UU<J)-U<J.IH0.5XSCDtJ> 

XX(J)*X(J,I H0 ,51SCA(  J  > 

YY<J)-Y<  J,I  )*0.5XSCB(J) 

CALL  PRIME (UU,UU, XX, YY, SCC, SCO, I,XFLAG) 

DO  14  J-l. 100 
SCC(J)-DRXSCCIJ) 

SCDU)-DRXSCD(J) 

R-R+0.5XDR 

RINU-1.0/R 

DO  15  J* 1,100 
UU<J)-U<J,I >+SCE<J> 

UU(J)-U<J,I )+SCF( J ) 

XX(J)-X(J,I)»SCC(J) 

YY(J)*V(J,1  )+SCD(  J ) 


CALL  PRinE<UU,UU.XX,W,SCE,SCF,I,IFLAG> 

DO  16  J-l, 100 

X<J.m>-X(J.I)4<SCG<Jl42.0*SCA(J)48.OXSCC(J>4DR*SCE(JmSIXINV 
YIJ,l*l  >-Y(J,I )4 (SCH( J)+2,0XSCB(J)t2.0XSCD(J)+DRXSCF(J))XSIXINU 

CONTINUE 


c  CALL  P«re<U<1.5),ga,S).X(l,5).Y(l,S),XP(l,fe  YP(1.S>.S,IFLAG> 
CXXSXXXSSXXXXXtXXtXXXXXXXXXXXXXXXXtXXXXXXXtttXXXXStXXXXXM 
c  "T  ™,|.£sTNrISI  fob  the 

DO  100  J-6.NDR 


R-R40R 
RZNU  -  1.0yR 

e 

•  DO  80  1-1,100 

X(I,J)-X(I,J-1 >4FACT* (55.0XXP(I , J-l )-59.0XXP( I, J-2) 

&  437.OJXP<X,J-3)-9.0»XP<I,J-4>> 

V(I,J)-Y< I, J-l )4FACTX(S5.0XYP(I,J-1 )-59.0XYP(I,J*8) 
t  437.0XVP(I,J-3>-9.O*VP(I,J-4>> 

U(I«J)*U(I, J-l  )4FACTX(55.0XX( I, J-l )-S9.0XX(I,J-2) 

1  437.0xX(I,J-3)-9.0XX(I,J-4>) 

WU,J)-W<I,J-l  )4F ACT*(55.0XVa, J-l J-59. 0*V<I,J-8> 

1  ♦37.0*Y(I,J-3)-9.0*Y<I,J-4)) 

80  CONTINUE 

C 

r  CALL  PRinE(U(l,J),U(l,Jl,X<l,J),Y(l,J),XP(l,J),VP(l,J),J,IFLAG> 

Cllttmtt  THE  CORRECTOR 
101  CONTINUE 

c 

DO  40  1-1,100 

SCA(l)-U(l.J-lHFACT<(9.0<X(I,J)»ig.0fX(I,J-l) 

t  -s.0xxu,j-e>+xu,j-3i> 

SCBII )-U(I, J-l )4FACTX(9.0XV( I, J )+19.0XV(I, J-I > 
l  -S.0*Y(I.J-2>4Y<I.J-3>) 

40  CONTINUE 
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DO  30  1*1,100 

Xd,J)*Xd,J-l  >*FACTS<9.0IXPd,J>*19.6SXPa.J-l> 
l  -S.0SXPd,J-2>+XPd.J-3>) 

y(i,j)*Ya,j-i)4FACT*(9.o*ypci,j)4i9.o*vp(i,j-i> 
t  -S.0aYP(I,J-2)4VP(I,J-3)) 

30  CONTINUE 

c 

CPU  PRINE(SCA,  SCI,  Xd,J>.YU,J>,XPU,J>.YPU,J>.J,IFlAG> 
C 

Clttimttt  ERROR  CHECK 
C 

ICORR-ICORRM 

IF ( ICORR.GT . 1000 )  60  TO  162 
DO  41  1*1,100 

SCCU  >*ABSUU(  I ,  J  >-SCAd  >  >/U(  X,J  >> 

41  SCD<I>*ABS((U(!,J)-SC»(ln/U(X,J>> 

DO  42  1*1,100 

Ud,  J  )*SCA( I ) 

42  U(I,J)*$CB(I> 

DO  43  1*1,100 

IF ( SCC (I >.LT. ERROR)  GO  TO  43 
GO  TO  101 

43  CONTINUE 

DO  44  1*1.100 


C 

100 

c 

102 

C 

CMS 

C 


603 

602 


600 


601 


C 

CMS 

c 


lF($CDd  >.LT. ERROR)  60  TO  44 

CO  TO  101 

CONTINUE 

CONTINUE 

CONTINUE 

USE  ANALYTIC  FIT  AT  SHEATH  BOUNDARY  SSSSSMSSSSSSSSSSMSSSSSSS 

DRSING-SNCLCDR) 

DO  601  1-1,100 

IF((DRSRS(1)).GE.(DRSRA<I)>)  GO  TO  603 
CO  TO  600 

IF(RA( I >.£0.0.0)  GO  TO  600 

UN(I)-1.0 

UN( I  )*1 .0 

XM( I  )*0.0 

YN(I )-DRSING/ABS< DRSING) 

UN.")1(I)*1.0 
UNN1C1  )*1.0 
XNMld  )*0.O 
VHni(I)*YN(I) 

RNNld  )-SNGL<RC> 

GO  TO  601 

IF(IFLAC(I ).LE.7>  GO  TO  602 
UNC I  )*Ud,  IFLACd  )-l) 

UNd  )*U(  I,  IFLAGd  )-l ) 

XN(  I  )*Xd,  IFLACd  )-l  ) 

VN< I  )*YCI,IFLAG( I )-l ) 

UNilld  )*Ud,  IFLAG! 1  )-2) 

L'Nflld  )*U(I, IFLACd  )-2) 

XNN1(  1  )*X(I, IFLACd  )-2) 

YNi11d)*Y(I, IFLACd  )-2) 

RNDld  )*SNGL<RC+DRSdFLAG(I)-2)) 

CONTINUE 

CALL  BOUND! DRSING, AC, BC,U0,5) 

SCALE  SELECTED  VARIABLES  INTO  CGS  UNITS  SSSSSSStSSSSSSMS 


USCA-3.0E2/EBYNC2 
RSCA • 1 . 0/ ( B0SEBYNC2 ) 

DO  210  1*1,100 
A< I  )»  AC(I)/ASCALE 
B(D*  BC(  I  )/BSCALE 
RS( I  )*RS< I )*RSCA 
RA(  I  )*RA(  I  XRSCA 
U (I , NDR )  *U (I , NDR ) SUSCA 
VPd,NDR)*Y(I,NDR)SBO 
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IF(LAZ)  CO  TO  210 

VP ( I , NDR )  •  YP  d  , NDR  )*Ud  , NOR )8B0/RANODE 
210  CONTINUE 
C 

cimtmmm  the  output  m««»mn»mmnmtmtmmmin 

c 

BOOTH  >  FXBOXO.BOl 
URITE (6,999)  BCATH 

999  FORMAT ( '  1 ' , /// , 20X , '  FOR  A  MAGNETIC  FIELD  NEAR  THE  CATHODE  Of', 

1  £18. S, '  KILOGAUSS  »',//> 

IF(LAZ)  GO  TO  291 
URITE (6, 1000) 

1000  F0RMAT<3X,'I',9X.'A',11X,'B',14X,  'RS' , 10X, 'RA'.llX, 'E(RA)', 

&  7X,'U<RA)',llX.'U<RA)',7X,'URA',//> 

DO  2e0  1*1.100 

200  URIT£(6. 1001 )  I.Ad  ),8d  ),RSd  >,RAd  ),Xd, NDR), Ud, NDR), U( I.NDR) 
1  .URA(I) 

1001  FORMAT!  IX.  I3,4(3X,2E12,4),/ > 

CO  TO  202 

201  IF ( LCIN )  CO  TO  203 
URITE(6. 1002 ) 

CO  204  1*1,100 

204  URITE(6.1O01 )  I.  A< I  ),B( I  ),RS(I  ),RA(  I  ),Xd,NDR),Ud,NDR) 

&  ,YP( I .NOR ),URA( I ) 

1002  FORMAT! 3X, ' I ' ,9X, ' A' , 1 IX, 'B ' , 14X, 'RS'.lOX, 'RA',11X,'E<RA)' 

&  ,7X, 'U(RA)', 10X, 'BT(RA)' ,6X, 'BTRA' ,// ) 

CO  TO  202 

203  URITE(6, 1003 ) 

DO  205  '1*1,100 

205  URITE(G, 1001 )  I,A(I),B(I  ),RS(I ),RA(I ),X( I ,NDR  ),U< I.NDR ) 

1003  FORMAT (2X, '  I ' . 12X, 'A' . 16X, 'B' , 15X, 'RS',13X, 'RA',12X, 
l  '£<RA)',8X,'U<RA)',//) 

202  CONTINUE 
C 

cmtltm  CALCULATION  OF  SAMPLE  B -FIELD  PROFILES 
C 

IF(LAZ)  GO  TO  251 
R  aRC 

DO  250  1*1, NOR 

R*R*DR 

RINU-l.O/R 

VPd,l)*(Vdl,I )+Udl,I )8RINV)  *  BO 
VP(2,I)*(V(I2,1)4U(I2,I)SRINU)  *  BO 

250  YP(3.I)-<Vd3.mud3,I)XRINV>  *  BO 
GO  TO  253 

251  CONTINUE 

DO  252  1*1 . NDR 
YP(l.I)  •  B0  *  Ydl.I ) 

YP<2, I )  -BO  *  Y(I2,I) 

252  VP(3,I)  •  BO  X  YCI3.I) 

253  CONTINUE 
C 

0X8  SCALING  OF  THE  POTENTIAL  X88X8888X8888SSS88XSX8888XXS8XS 

C 

00  260  I« I.NDR 
U(Il.I)  •  Udl.I)  X  USCA 
U( 12, I )  ■  U< 12, I )  X  USCA 
260  U(I3.I)  •  U( 13,1 )  X  USCA 
C 

URITE16, 1006)  11,12,13 

1006  FORMATdX, /✓///, IX, 'PHI  AND  B -FIELD  PROFILES  FOR  (A,B)  CHOICES', 
1  14, ','.14,',  AND'. 14,' 

DO  300  I* I.NDR, 2 

300  URIT£(6,1005)  I.Udl ,1  ),YPt  1.1  ),U(I2,I  ),YP(2.I ),Ud3,X),YP(3,I ) 
1005  FORMATdX,  13,31 3X.2E14.5)) 

CXXX  COMPARE  VALUES  AT  ELECTRON  SHEATH  BOUNDARY  XXXXXXXXXXXXXXXXX 
C 

CALL  HATCH(RSTAR,URA, YSTAR.XSTAR.BC ) 

C 

URITE(6, 1010 ) 


NRL  MEMORANDUM  REPORT  4654 


let*  F0RHAT(1X,/////,1X,  'HATCHING  VALUES  AT  SHEATH  BOUNDARY',//, 
4  3X,'I',9X,'RS',13X,'US1',8X. 'US2',12X.'ES1\ 

4  8X.'ES2',12X,'US1',0X,'YS1',9X,'US2'.//) 

cm  SCALE  RSTAfi  INTO  CENT1NETERS  tXtUtttttWtUtttttttSSSt 

DO  650  1-1,100 
650  RS( I  )-RSTAR( I  )SRSCA 
C 

DO  700  1-1,100 

70*  URITE(6, 10041  I,RS(1 ),USTAR(1 ),US( I  ),XSTAR(I  ),ES( I ), 

4  USTAR(  1 ),  YSTARC 1  >,US(  1 1 

c 

1004  FORMAT! IX,  I3,2X,E12.4,2(3X,2E12.4),3X,3£12.4,/) 

IF(IT.LT.NF)  GO  TO  400 

c 

STOP 

C 

END 


C 

C 

C 

c  subroutine  bound 

c - 

c 

SUBROUTINE  BOUND(DR.A,B,U0,XTER) 

LOGICAL  LAZ.LC1N, LIONS 
C 

DIMENSION  A<100),B(100),R<100),SCR1(100> 

COfinON  /FIT/  UN(100),UN( 100 ) ,XN( 1C0 ) , YH< 100 ) ,UNN1 ( 100 ),UNN1 ( 100), 

4  XNN1  ( 100 ) ,  YNP.l  ( 100  J.RNfU  ( 100  >,R$TAR(  100  ),USTAR( 100  ) ,USTAR(  100  >, 

4  XSTAR< 100 ),YSTAR< 100 ) 

COrtflON  /SWITCH/  LAZ.LCIN.  LIONS 

COMMON  /ITCHY/  A1 ( 100), B1 ( 100), A2( 100 ),B2C 100 > 

C 

EPS-1.0 

IF(LAZ)  EPS-0.0 
F  -  1.0 

IF ( . NOT.LCIN)  F— 1.0 
C 

cm  hake  zeroeth-order  guesses  mmmmmmmtmmmm 
c 

DO  1  1-1,100 

R(I)-DR 

USTAR(I)-UN<I) 

USTARd)-UN(I) 

XSTARt I  )-XN( I ) 

i  ystarcd-ynu) 

c 

Cm  ITERATE  FOR  SOLUTIONS  tt 

DO  2  IT- 1, ITER 
C 

CM*  CALCULATE  A1  AND  11  Mf 
C 


3 


C 


DO  3  1-1,100 

A2< I )-F>4.0tA(I )/(3.0<(R( 1  )+RNNl( I  ) ) ) 

B2( I  )-F((2.SUSTAR< I )*YSTAR< I )-2.IXSTAR<X )*(1 .+USTARU ))) 


20 

C 


DO  20  1-1.100 
IF(B2(I).GT.0.0)  GO  TO  2* 
B2( 1 1-1.0E20 
CONTINUE 


DO  4  1-1,10* 

Al( 1 )•  A2< i )*< 1 ,04U$TAR( X ) )/S0RT(B2(X ) ) 


ooooooo 
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4  ASd)*USTARd)/'SCRT<B2d>> 

C 

C*M  CHECK  FOR  NEGATIVE  VALUES  OF  'USTAR'  tumtttttutilttm 
C 

DO  30  1-1,100 

IF ( USTAR ( I ), GT. 0.0. AND, USTAR ( X ),LT,U0)  GO  TO  30 
USTARC I  )•  (  1 .0E-3XU0 
30  CONTINUE 
C 

Ct«  CALCULATE  R-CORRECTED  *** 

C 

DO  5  1-1,100 

AS (I)-  SORT  <  UST AS?  d  )  *  <  2 . O+USTAR (I))) 

S  B2(  I  )•  1.0  ♦  USTARd) 

C 

DO  6  1-1,100 

6  RSTARd  )•  (A2d)S(1.0HINHld))-UNnid)SBSd))  / 

1  (VSTAR(  I  )*B2( I )  -  XSTARd  )*A2d>) 

C 

DO  14  1-1,100 

14  IF(<RSTARd)*F>.GE.O.O>  R(I)-RSTARd) 

C 

CtM  CALCULATE  AS  AND  B2  Ht 
C 

DO  7  1-1,100 

7  SCRKI)—  0.5/<RdHRNHld)> 

C 

DO  8  1-1,100 

A2d  l-SCRl ( I  )>(XSTAR( I )*8( I )/SORT(UO-USTAR(I ) ) > 

8  B2(  I  )-SCRl  ( I  )t(VSTAR(  I  )*2. 0*EPS*USTAR( I  XSCRl (I  > ) 

C 

CM*  CORRECT  THE  PREDICTED  UALUES  *«* 

C 

DO  9  1-1,100 

XSTARd  )-XNMl( I  )»1 ,5*A1  ( I  )*SORT(F*R(I )  >+2.0*A2d  )8R(I ) 

9  YSTARt  I  )-YNf11  ( I  )*1 .5SB1  ( I  )*SORT (F8RCI  >  >+2.0*B2( I  )*R(I  ) 

C 

DO  10  1-1.100 

USTAR i  I ) -UNI1 1  d  HR <1  )* ( XSTARd  )-Al  d  )*SORT ( F*R(I )  >-A2(  I  >SR< I )  ) 

10  USTARt I )-UNHl ( I  )*R< I )t< VSTAR< I )-Bl 1 1  )*SORT(F*R( I ) >-B2t I )*R( I ) ) 
C 

2  CONTINUE 
C 

DO  12  1-1,100 
12  RSTAR<  1  )-RNNld  )+Rd  ) 

C 

RETURN 

C 

END 


SUBROUTINE  PRIME 


SUBROUTINE  PRIME<U,U,X,V,XP,VP,IT,XFLAO> 

DIMENSION  S(lOO),U<10O),U(lOO).X(lOO),Y(lOO),XPaOO),YP(lOO) 

COMMON  /ITCHV/SCA(  100 ),SCB(  10O),SCC(  100>,SCD(  100)  ' 

COMMON  /COEF/  A(100),B(100),RS(100),RA(100),AFLAGdOO)» 

1  R,DR.RC,RAN0DE,RINV,U0,B0,EBYMC2,EMRAT,C,GAPLEN 
COMMON  /SUITCH/  LA2.LCIN, LIONS 

LOGICAL  LA2.LCIN, LIONS 
INTEGER  !pLAG(100) 

REALS8  R,DR,RC,RANODE 
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C 

c ttttttxtlt  FIND  THE  ELECTRON  SHEATH  RADIUS  XXXXXXXXXttXXXXXXXXXXXXXXX 
C 

DO  1  I -1.180 

1  s(n-2.0tu(i)4rum«u<ins(u(i)-u(i>) 

DO  2  1-1,100 

IF (SI  I  ).GT.l .8E-8 )  GO  TO  Z 

SCI )-1.0 

ACD-0.0 

1F( IFLAGIX ).EO.0)  XFLAGtl )-XT 

2  CONTINUE 
C 

DO  7  1-1,100 
7  SCBf I )-A( I )tRS(I ) 

C 

DO  3  1-1,100 

1F(SCB(1 1.NE.0.0)  GO  TO  3 
RS(I  )-R 

3  CONTINUE 
C 

cmmtm  find  the  'anode  radius'  ttiuimmmmmmummu 

DO  4  1-1,100 

IF(U(I).LT.U0)  GO  TO  4 

A1D-0.0 

KD-0.0 

AFLAGC1 )-0,0 

X<I)-0.0 

Y< 1  )-0.0 

4  CONTINUE 
C 

DO  S  1-1,100 

5  SCA( I  )-AFLAG( I  )-RA( Z ) 

C  • 

DO  6  1-1,100 

IF(SCA( I1.NE.0.0)  GO  TO  « 

RA(I)-R 

6  CONTINUE 
C 

DO  10  1-1.100 

10  SCA(l)-  1.0/SQRT(S(D) 

DO  11  1-1.100 

11  SCB(I)-A(I)SSCh(I) 

DO  12  1-1,100 
XP(I)-SCB<I)*(1.0HJ(I)> 

12  VP(  I )«SCB( I )*U( I ) 

DO  13  1-1,100 

13  SCA(I)  •  U0-U(I) 

DO  14  1-1.100 

IF  <SCA( 1 ).GT.0.O)  GO  TO  14 
SCA(I)-1.0 

14  CONTINUE 

DO  IS  1-1,100 

15  SCC(I>-1.0'SORT<SCA(I>) 

DO  16  1-1,100 

16  SCA(I)-XPU>-B(I)*SCCm 
DO  17  1-1,100 

17  SCB(I )  •  YP< I ) 

IF(LAZ)  GO  TO  19 
DO  18  1-1,100 

18  SCBU)-SCB<I>+AFIA6(I)SU<X)SRXNU 

19  CONTINUE 

DO  20  1-1,100 
SCC( I )-SCA( I )-X( 1 ) 

20  SCD( I )-SCB( I )-Y(I ) 

DO  21  1-1,100 
XP(I )*SCC(I ttRINU 

21  YP( I )*SCD( X HRINU 
C 

c 


RETURN 


QOOOOOO 
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SUBROUTINE  NATCH 


C 

C 


C 

cm 

c 


1 

2 


3 

C 

cttt 

c 


7 

6 


4 

9 

C 

8 

C 

c*« 

c 


SUBROUTINE  NATCH ( RSTAR, URA, VSTAR, XSTAR, BC > 

LOCICAL  LAZ.LCIN, LIONS 
REALM  R, DR. RC.R ANODE 

CONNON  /COEF/  A<  100). 8(100  >,RS<  100), RAUO0),AFLA6<100># 
l  R,DR,RC,RANODE,RINV,U0,B0 

COIWON  /SWITCH/  LAZ, LCIN, LIONS 

CONNON  /ITCHV/  SCAC 100 ),SC8< 100 >,SCC< 100 >,SCD(100) 

D INENSI ON  RSTAR (1 00 ) . UR A ( 1 00 ) . VSTAR ( 1 00 ) , XSTAR ( 100 ) , 
l  BC< 100  >,RAT< 100 >,RATINV< 100 ),G( 100  ),BETA( 100 ), 

0  BETAPC 100) 

be  tern  ine  needed  constants  immmmmmmmm 

RASING*  SNGL  C  RANO  DE ) 

DO  1  1*1,100 

RAT( I ) • RAS I NG/RSTAR ( 1 1 

C(I)*1.0 

DO  2  1*1,100 

RATINUd )*1 ,0/RAT( 1 ) 

DO  3  1*1,100 

IF (RATINV< I ) .GT.0.0)  G< I >*ALOG(RATINUd ) ) 

CONTINUE 

NATCH  THE  HAGNETIC  VARIABLE  tttttttt tttttttttttttttttttt 
IF (LAZ)  GO  TO  4 
DO  S  1*1.100 

SCA( I  )*  2.0<URA(I )  -  YSTARd  )<<RASING-RSTAR< I )<RATINV( I ) > 

DO  6  1*1.100 

IF(RAT( I  ).NE. 1 .0 }  GO  TO  7 
GO  TO  6 

SCB(I)  ■  SCA< I )/  <RAT( I )«RATINV( I ) ) 

CONTINUE 
GO  TO  8 

CONTINUE 
DO  9  1*1.100 

SCB(I)  •  URA<D<RAT<  I  >/B0 
CONTINUE 

NATCH  THE  ELECTRICAL  VARIABLES  »«X»»»M»«»»»t»»K 
IF(. NOT. LIONS)  GO  TO  20 


10 


00  10  1*1,100 

BETA(I)  •  G( I ><(  1.0  -  G( I )<(  0.400  -  C(I)t(  9.16S6E-2 
4  -  G< I  )<(  1.4242E-2  -  G< I )<(  1.6793E-3  -  G(I><(  1.6122E-4 

1  -  C(I)t(  1.2935E-5  -  G(I)<(  8.87G9E-7  -  G<I)t<  S.4619E-8 

0  -  G(  I  )<(  2.9484E-9  -  G(I >*(  1.3603E-19  -  G(I)<<  7.U01E-13 

4  -  GC I >*(  2.6644E-13  ♦  C<I)<  1 .2S26E-1S) ))))))))))) ) 

CONTINUE 


11 


00  11  1*1,100 

BETAP(I)  •  (  1.0  -  C(I)«  0.8  -G<I)«  2.7SE-1  -  G(D« 

4  S.S9G8E-2  -  G<1)<(  8.3964E-3  -  C(l)t(  9.6733E-4 

4  *  G(  I  )<<  9.0544E-5  -  G(I»(  7.101SE-6  -  G(I)«  4.91S7E-7 

4  -  G( I  )<(  2.9484E-8  -  G(I)«  1.49G3E-9  -  G(D<(  8.S321E-12 

4  -  G< I )<(  3.4637E-12  ♦  G(I><  1 .7S36E-14 )))))))))))) ) /RSTAR (I ) 

CONTINUE 


30 
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DO  12  Is!  100 

12  SCC<I>  •  Ue  -  <ABS<8.6*RSTAR<n*BETA<I>*BETA<I)*BCU>>/4.6> 

fc  tS(2.e/3.0> 

DO  13  1*1  100 

13  SCD(1>«  -<2.0/3.e)t(U0  -  SCC(I))m.e/RSTAR(I) 

t  *  2.0*BETAP(I)/IETA(I)) 

CO  TO  H 

CZ$  CONTINUE 
C 

DO  15  I-t.160 

15  If<RAT(I  I.GT.0.0)  C(I)*ALOG<RAT(I>> 

DO  16  1«1,1M 

16  5CC(I)  -  DA  -  XSTARt 1  )*RSTAR< I )*C( I ) 

C 

14  RETURN 
C 

END 


Appendix  B 

SAMPLE  SOLUTION  FOR  A  RADIAL  DIODE 


In  order  to  provide  an  example  of  the  iterative  "shooting”  technique  solution  process,  the  radial 
diode  depicted  in  Fig.  6  was  modeled.  Two  megavolt  operation  of  such  a  geometry  should  make  the 
case  interesting  to  those  in  the  light  ion  beam  research  area.  The  desired  configuration  was  specified  by 
the  following  choice  of  variables: 

a.  LAZ  -  FALSE., 

b.  LCIN  -  TRUE., 

c.  LIONS  -  TRUE., 

d.  RC  -  5.0, 

e.  RANODE  -  5.5, 

f.  U0  *■  2.0E6,  and 

g.  B0  -  2.0E4. 


As  for  the  numerical  factors,  NDR  was  left  at  500  and  ERROR  was  chosen  as  0.02. 

'  3  2cm——]  R] 

1.0 1— 


Fig.  6  —  The  radial  diode  tesl  case 
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A  quasi-static  computer  simulation  of  this  diode  has  been  completed25  and  the  numerical  results 
indicate  an  electron  current  density  at  the  cathode  of  approximately  3.8  x  104  amps/cm2  with  a  comple¬ 
mentary  ion  current  density  of  about  1.3  x  104  amps/cm2.  Due  to  the  geometric  electric  held  enhance¬ 
ment  at  the  anode  surface,  it  can  be  assumed  that  this  ion  current  is  too  high.  In  the  simulation,  the 
axial  magnetic  held  strength  at  the  cathode  surface  was  found  to  go  slightly  negative.  Such  held  rever¬ 
sal  cannot  be  expected  if  the  electron  and  ion  currents  are  signihcantly  reduced.  Given  the  high  vol¬ 
tages  of  this  test  case,  the  approximate  solutions  presented  in  Figs.  3  and  4  do  not  provide  any  addi¬ 
tional  assistance  in  choosing  reasonable  guesses.  As  a  hrst  test  for  RADBER,  therefore,  it  was  decided 
to  search  the  following  parameter  ranges: 

(a)  (flz)cathode  “  4,  5,  6,  7,  and  8  kilogauss. 

(b)  Ten  equispaced  values  of  2.0  x  104  <  j  <  4.25  x  104  amps/cm2. 

(c)  Ten  equispaced  values  of  0.3  x  104  <  J  <  0.75  x  104  amps/cm2. 

In  examining  the  numerical  results  from  the  program  it  is  important  to  note  that  in  cases  such  as  this 
where  the  electron  sheath  boundary  will  lie  close  to  the  anode  surface,  the  matching  of  values  at  p  - 
p, 4  is  extremely  unreliable  due  to  the  sharp  discontinuities  in  U"  and  W"  at  p  -  ps.  As  a  general  rule, 
therefore,  one  must  perform  the  matching  tests  at  p  -  ps-  After  good  matches  are  found  there,  the 

corresponding  tables  for  p  -  p  A  may  or  may  not  add  corroboration  but  they  can  never  be  used  to  over¬ 

rule  the  p s  test  results.  Similarly,  the  sample  <t>  and  Bz  profiles  are  to  be  trusted  completely  only  inside 
the  electron  sheath. 

An  examination  of  the  matching  results  for  this  first  parameter  search  yielded,  as  might  be 
expected,  only  a  few  promising  candidates.  By  far  most  of  the  values  to  be  compared  at  p  -  p  5  differed 
by  at  least  factors  of  two.  From  this  mass  of  data,  however,  the  best  matches  were  extracted  for 
presentation  in  Table  2.  The  closest  fit  seems  to  be  reached  for  -  7  kG,  j  -  2.75  x  104  and 

J  -  0.35  x  104.  The  next  logical  search  region  seems  to  be: 

33 
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TABLE  2.  SAMPLE  RADIAL  DIODE  OUTPUT 


BZ 

JE 

JI 

(US 'IN 

(US 'OUT 

(XS )  IN 

(XS)OUT 

(US)IN 

(US)OUT 

5.0 

2.50 

0.40 

2.36 

3.89 

0.552 

0.559 

3.21 

5.54 

5.0 

4.00 

0.55 

2.95 

3.89 

0.443 

0.688 

3.82 

5.52 

6.0 

2.70 

0.40 

2.85 

3.76 

0.703 

0.956 

3.72 

5.25 

6.0 

3.50 

0.45 

3.36 

3.91 

0.532 

0.340 

4.25 

5.58 

7.0 

2.25 

0.30 

3.23 

3.87 

0.855 

0.612 

4.11 

5.44 

7.0 

2.75 

0.35 

3.48 

3.86 

0.772 

0.694 

4.36 

5.41 

7.0 

3.50 

0.40 

3.74 

3.85 

0.293 

0.776 

4.63 

5.33 

8.0 

2.50 

0.30 

3.66 

3.73 

0.906 

0.871 

4.55 

5.08 

UHERE  IN  THIS  TABLE  AS  UELL  AS  IN  TABLES  3,  4,  AND  5: 

BZ  IS  IN  KILOGAUSS, 

JE  AND  JI  ARE  IN  UNITS  OF  10,000  AMP'CMSG,  AND 

US,  XS,  AND  US  ARE  ALL  DIMENSIONLESS. 

(a)  (Bz) ctlh  -  6.5,  7.0,  and  7.5  kG. 

(b)  Ten  equispaced  values  for  2.0  x  104  <  j  <  3.8  x  104. 

(c)  Ten  equispaced  values  for  0.25  x  104  <  J  ^  0.43  x  104. 

The  RADBER-generated  matching  values  at  p  -  ps  provided  the  results  listed  in  Table  3.  The  best 
overall  answers  seem  to  be  for  (fl,)CIth  “  7.5  kG,  j  ■*  2.80  x  104,  and  J  -  0.33  x  104.  For  those 
parameters,  the  worst  match  occurs  for  Ws-  However,  the  jump  in  those  values  of  Ws  compared  to 
their  neighbors  in  the  table  suggests  that  the  program  is  searching  in  a  sensitive  region  of  parameter 
space.  Such  "sensitivity"  often  implies  that  one  is  close  to  the  correct  answer  in  the  search.  To  further 
delineate  the  correct  range,  it  seems  reasonable  to  check  adjacent  (fl2)CIIh  values.  For  the  third  test 
search,  therefore,  the  following  ranges  are  chosen: 


(a)  “  7.25,  7.50,  7.75,  and  8.00  kG. 

(b)  Ten  equispaced  values  of  1.7  x  104  <  j  <  3.5  x  104. 


(c)  Ten  equispaced  values  of  0.25  x  104  <  J  <  0.43  x  104. 
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TABLE  3.  SAMPLE  RADIAL  DIODE  OUTPUT 


BZ 

JE 

JI 

(US  JIN 

(US j OUT 

t  >:s  jin 

iXSJOUT 

(US  JIN 

(USJ0UT 

6.5 

2.00 

0.29 

2.98 

3.89 

0.790 

0.526 

3.85 

5.51 

6.5 

2.20 

0.31 

2.81 

3.86 

0.788 

0 . 635 

3.68 

5.42 

6.5 

3.00 

0.39 

3.38 

3.89 

0.688 

0.604 

4.26 

5.50 

6.5 

3.20 

0.41 

3.43 

3.86 

0.643 

0.732 

4.32 

5.43 

7.6 

2.00 

0.29 

2.76 

3.60 

0.350 

0.981 

3.63 

4.95 

7.0 

2.20 

0.31 

2.94 

3.66 

0.360 

0.962 

3.81 

5.03 

7.0 

2.40 

0.33 

3.09 

3.69 

0.357 

0.959 

3.97 

5.08 

7.0 

2.60 

0.35 

3.21 

3.70 

0.341 

0.975 

4.09 

5.09 

7.0 

2.80 

0.35 

3.55 

3.39 

0.732 

0.559 

4.44 

5.50 

7.0 

3.00 

0.37 

3.60 

3.86 

0.674 

0.698 

4.49 

5.42 

7.5 

2.00 

0.27 

3.06 

3.67 

0.938 

0.390 

3.93 

5.01 

7.5 

2.20 

0.29 

3.23 

3.71 

0.934 

0.384 

4.11 

5.07 

7.5 

2.40 

0.31 

3.36 

3.72 

0.914 

0.399 

4.24 

5.09 

7.5 

2.60 

0.33 

3.46 

3.71 

O.380 

0.934 

4.35 

5.03 

7.5 

2.80 

0.33 

3.75 

3.87 

0.667 

0.641 

4.65 

5.43 

7.5 

2.80 

0.35 

3.54 

3.69 

0.335 

0.938 

4.43 

5.03 

7.5 

3.00 

0.35 

3. 78 

3.82 

0.530 

0.784 

4.67 

5.31 

The  best  matches  are  listed  in  Table  4.  Clearly  the  most  tantalizing  results  are  those  for  the  8.0  kilo- 
gauss  case.  For  that  fixed  value  of  (flz)calh  a  dose  match  for  each  of  the  three  variables  tested  can  be 
found  within  the  given  (J.J )  matrix.  More  importantly,  all  of  the  matches  fall  on  the  line  j  -  2.10  + 
10  (J  -  0.27)  -  107-  0.60.  Two  possibilities  suggest  themselves;  either  the  solution  falls  on  a  similar 
line  for  an  adjacent  magnetic  field  value  or  it  lies  on  an  adjacent  (J.J) -line  unresolved  by  the  coarse 
mesh  but  at  the  same  field  value.  The  former  alternative  will  be  tested  first  by  examining  the  same 
(J.J)  mesh  for  the  values,  (fl2)c„h  -  7.9,  8.1,  8.2,  8.3,  8.4,  8.5.  The  major  results  are  listed  in  Table  5. 
The  coalescence  of  circled  solutions  indicates  a  finer  scale  search  in  the  region  delineated  by: 

(a)  («,)„. h  -  8.3,  8.4,  8.5,  8.6,  8.7  kG. 

(b)  Ten  equispaced  values  of  1.65  x  104^  j  <  2.55  x  104. 

(c)  Ten  equispaced  values  of  0.225  x  104  <  J  <  0.315  x  104 

No  tabulated  results  are  necessary  for  this  run  since  a  "hit"  was  scored  in  the  p  -  ps  matching  results 
for  the  parameter  combination  (fl,),,,!,-  8.7  kG,  j  -  2.25  x  104  amps/cm2,  and  J  -  0.275  x  104 
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TABLE  4 .  SAMPLE  RADIAL  DIODE  OUTPUT 


BZ 

JE 

JI 

(US)IN 

(US)OUT 

(XS)IN 

(XS)OUT 

(US)IN 

(US)OUT 

7.25 

2.10 

0.29 

3.00 

3.66 

0.899 

0.928 

3.87 

5.02 

7.25 

2.30 

0.31 

3.16 

3.70 

0.895 

0.924 

4.04 

5.07 

7.25 

2.50 

0.33 

3.29 

3.71 

0.878 

0.939 

4.17 

5.09 

7.25 

2.90 

0.35 

3.68 

3.86 

0.678 

0.671 

4.57 

5.42 

7.25 

3.10 

0.35 

3.81 

3.91 

0.321 

0.338 

4.70 

5.58 

7.25 

3.10 

0.37 

3.70 

3.82 

0.605 

0.807 

4.60 

5.32 

7.50 

2.50 

0.31 

3.56 

3.83 

0.859 

0.723 

4.45 

5.33 

7.50 

2.70 

0.33 

3.63 

3.81 

0.801 

0.792 

4.52 

5.27 

7.50 

2.90 

0.33 

3.83 

3.89 

0.459 

0.501 

4.73 

5.52 

7.50 

2.90 

0.35 

3.68 

3.77 

0.734 

0.879 

4.57 

5.19 

7.50 

3.10 

0.35 

3.83 

3.85 

0.366 

0.726 

4.72 

5.37 

7.50 

3.10 

0.37 

3.71 

3.73 

0.662 

0.975 

4.60 

5.09 

7.50 

3.30 

0.37 

3.82 

3.78 

0.288 

0.888 

4.71 

5.22 

7.50 

2.90 

0.39 

3.22 

3.26 

0.908 

1.374 

4.09 

4.40 

7.50 

3.10 

0.39 

3.53 

3.52 

0.813 

1.207 

4.41 

4.75 

7.75 

2.10 

0.27 

3.29 

3.72 

0.972 

0.841 

4.17 

5.07 

7.75 

2.30 

0.29 

3.43 

3.73 

0.950 

0.855 

4.32 

5.10 

7.75 

2.50 

0.31 

3.54 

3.73 

0.910 

0.891 

4.43 

5.08 

7.75 

2.70 

0.31 

3.83 

3.87 

0.627 

0.610 

4.72 

5.43 

7.75 

2.70 

0.33 

3.62 

3.70 

0.856 

0.947 

4.51 

5.04 

7.75 

2.90 

0.33 

3.84 

3.82 

0.517 

0.762 

4.74 

5.30 

8.00 

2.10 

0.27 

3.26 

3.59 

1.011 

0.958 

4.14 

4.82 

8.00 

2.30 

0.29 

3.40 

3.61 

0.994 

0.977 

4.29 

4.85 

8.00 

2.50 

0.31 

3.52 

3.60 

0.960 

1.012 

4.41 

4.84 

8.00 

2.70 

0.33 

3.61 

3.58 

0.911 

1.063 

4.50 

4.80 

8.00 

2.90 

0.35 

3.67 

3.54 

0.849 

1.125 

4.56 

4.74 

8.00 

3.10 

0.37 

3.72 

3.49 

0.778 

1.197 

4.61 

4.65 

8.00 

3.30 

0.39 

3.75 

3.42 

0.701 

1.276 

4.64 

4.55 
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TAILS  S.  SAMPLE  RADIAL  DIODE  OUTPUT 


•z 

J€ 

J! 

(US  JIN 

(US  )0UT 

(XS  >IN 

(XS  >0UT 

(US  JIN 

(US  >GUT 

7.9 

2.39 

9.29 

3.42 

3.66 

9.977 

0.932 

4.39 

4.95 

7.9 

2.5# 

9.31 

3.53 

3.65 

9.949 

9.968 

4.42 

4.94 

7.9 

2.7# 

9.33 

3.61 

3.63 

9.889 

1.029 

4.59 

4.99 

7.9 

2.9# 

9.35 

3.68 

3.S9 

9.826 

1.08S 

4.57 

4.83 

7.9 

3.1# 

9.37 

3.72 

3.54 

8.754 

1.169 

4.61 

4.74 

7.9 

3.3# 

9.39 

3.7S 

3.47 

0.677 

1.241 

4.64 

4.64 

7.9 

3.5# 

9.41 

3.76 

3.49 

9.599 

1.133 

4.65 

4.52 

S.l 

2.1# 

9.27 

3.24 

3.53 

1.926 

9.998 

4.12 

4.71 

S.l 

2.3# 

9.29 

3.39 

3.55 

1.911 

1.918 

4.28 

4.75 

1.1 

2.5# 

9.31 

3.51 

3.55 

0.989 

1.953 

4.49 

4.74 

a.i 

2.7# 

9.33 

3.69 

3.53 

0.932 

1.192 

4.49 

4.71 

8.1 

2.9# 

9.35 

3.67 

3.49 

0.873 

1.163 

4.56 

4.65 

a.i 

3.1# 

9.37 

3.72 

3.44 

0.892 

1.233 

4.61 

4.56 

8.2 

2.1# 

9.27 

3.22 

3.47 

1.040 

1.035 

4.10 

4.61 

8.2 

2.3# 

9.29 

3.38 

3.59 

1.028 

1.956 

4.26 

4.65 

8.2 

2.5# 

9.31 

3.59 

3.59 

0.999 

1.991 

4.39 

4.64 

8.2 

2.5# 

9.33 

3.29 

3.16 

1.039 

1.314 

4.98 

4.18 

8.2 

2.7# 

9.33 

3.59 

3.48 

0.954 

1.139 

4.48 

4.61 

8.2 

2.9# 

9.35 

3.66 

3.44 

0.397 

1.198 

4.56 

4.55 

8.3 

1.7# 

9.2S 

3.91 

3.36 

1.045 

1.958 

3. 88 

4.44 

8.3 

2.1# 

9.27 

3.21 

3.41 

1.054 

1.069 

4.03 

4.51 

8.3 

2.39 

9.29 

3.37 

3.44 

1.045 

1.092 

4.2S 

4.55 

8.3 

2.59 

9.31 

3.49 

3.44 

1.018 

1.127 

4.38 

4.55 

8.3 

2.79 

9.33 

3.59 

3.42 

0.975 

1.174 

4.48 

4.52 

8.3 

2.39 

9.35 

3.66 

3.39 

0.918 

1.232 

4.55 

4.46 

8.4 

1.39 

9.25 

2.99 

3.29 

1.057 

1.088 

3. 86 

4.34 

8.4 

2.19 

9.27 

3.19 

3.35 

1.069 

1.101 

4.07 

4.41 

8.4 

2.39 

9.29 

3.35 

3.38 

1.062 

1.125 

4.23 

4.45 

8.4 

2.59 

9.31 

3.48 

3.39 

1.937 

1.169 

4.37 

4.45 

8.4 

2.79 

9.33 

3.58 

3.37 

0.996 

1.207 

4.47 

4.42 

8.5 

1.39 

9.25 

2.97 

3.23 

1.067 

1.117 

3.84 

4.24 

8.S 

2.19 

9.27 

3.17 

3.29 

1.982 

1.131 

4.05 

4.31 

8.5 

2.39 

9.29 

3.34 

3.32 

1.977 

1.156 

4.22 

4.35 

8.5 

2.59 

9.29 

3.99 

2.92 

1.069 

1.365 

3.88 

3.E8 

8.5 

2.59 

9.31 

3.47 

3.33 

1.056 

1.192 

4.35 

4.35 

amps/cm2.  For  that  choice,  RADBER  found  the  following:  rs  —  5.445,  (l/s)m  =  3-438,  (f/s)0Ui“ 
3.381,  (*s)in-  1.104,  Us)ou(  -  1.096,  OFs)in-  4.324,  and  ( M^s)oul  —  4.397.  It  is  interesting  to  note 
that  the  matching  at  p  -  pA  yielded  rs  *  5.445,  rA  -  5.497,  EA  -  -0.0053,  <t>  A  ■»  2.0  MV,  and  - 
5.435.  Only  the  value  for  WA  is  significantly  displaced  from  its  proper  value  of  5.600.  A  final  fine- 
tuning  run  is  called  for  not  so  much  to  obtain  better  values  for  the  parameter  guesses  but  rather  to 
obtain  useable  profiles  for  the  quantities  <Mr)  and  BAr).  For  this  final  run,  it  was  decided  to  pick: 

(a)  (B,)calh  -  8.675,  8.700,  and  8.725  kG. 

(b)  Ten  equispaced  values  of  2.05  x  104  ^  j  <  2.50  x  104. 

(c)  Ten  equispaced  values  of  0.255  x  104  <  J  <  0.300  x  104. 
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The  results  of  that  run  clearly  indicate  a  best  choice  of  parameters  as  l(£z)c4th>  -/.•/)!  ”  18.725  kG,  22.0 
kA/cm2,  2.75  kA/cm2!.  For  that  combination,  matching  at  p  -  ps  yielded  (L/S)m  -  3.306,  ({/s)oul  - 
3.268,  (Es) iB-  1.116,  (£s)out-  1.152,  OFs)in-  4.188,  OFs)out-  4.235,  and  rs  -  5.436  cm.  The  per¬ 
centage  errors  in  each  of  the  respective  quantities  are  -  1.2%,  A£  -  3.2%,  and  -  1.1%.  Given 
the  2%  error  margin  chosen  for  the  predictor-corrector,  it  is  mathematically  pointless  to  seek  a  closer 
match.  The  radial  profiles  of  the  electrostatic  potential  and  axial  magnetic  field  values  supplied  by 
RADBER  for  that  chosen  solution  are  plotted  in  Fig.  7.  Note  that  the  nature  of  the  B-field  profile  is 
qualitatively  identical  to  those  plotted  by  Ron,  et  a!.,  (see  Ref.  11).  On  the  other  hand,  the  E-field 
inferred  from  the  curve  of  in  Fig.  7  sharply  differs  from  that  cited  in  the  above  journal  article. 
The  difference  exists  because  the  presence  of  ions  was  not  treated  in  Ref.  11.  Ion  emission  at  the 
anode  surface  requires  that  the  electric  field  strength  go  to  zero  there. 


Fig.  7  —  Profiles  of  <4  and  B.  in  ibe  A-K  gap  of  ihe  radial  diode. 


These  RADBER  results  represent  a  very  significant  departure  from  the  results  of  the  full 
DIODE2D  computer  simulation  of  the  radial  diode  under  study.  Instead  of  a  slightly  negative  Bz  at  the 
cathode  surface,  RADBER  predicts  8.725  kG.  Instead  of  an  electron  current  density  of  38  kA/cm2, 
RADBER  finds  22  kA/cm2.  Finally,  and  most  dramatically,  we  are  left  with  only  2.75  kA/cm2  of  ion 
current  density  compared  to  DIODE2D’s  13.0  kA/cm2.  These  observed  differences,  bovver,  may  be 


attributed  to  two  mqjor  causes: 
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(1)  The  electric  field  near  the  anode  of  the  actual  simulated  radial  diode  is  enhanced  due  to  2-D  dis¬ 
tortions  which  cannot  be  modeled  by  RADBER’s  one-dimensional  treatment. 

(2)  The  D10DE2D  simulation  accurately  treated  the  azimuthal  magnetic  insulation  present  in  the 
actual  diode  as  well  as  the  axial  one.  Given  the  currents  of  over  a  megampere  flowing  through 
the  cathode  shank,  very  large  values  of  Be  could  be  expected.  This  additional  field  component 
could  easily  have  modified  the  results  of  the  simple  1-D  model  by  the  degree  shown.  As  an  addi¬ 
tional  reference  point,  simple  Child-Langmuir  bi-polar  analysis  predicts  an  electron  current  density 
of  just  under  50  kA/cm2  and  a  corresponding  J  of  1.2  kA/cm2.  It  is  interesting  that  both  densities 
are  off  by  factors  of  two  from  the  RADBER  results  but  in  opposite  directions. 

In  conclusion,  therefore,  RADBER,  has  given  a  reliable  order-of-magnitude  estimate  for  the 
current  flow  through  the  system.  It  also  gives  an  investigator  some  appreciation  for  the  degree  to  which 
combined  axial  and  azimuthal  insulation  can  modify  a  simple  axial  insulation  prediction.  Of  course, 
RADBER  can  also  be  used  to  predict  a  purely  azimuthal  insulation  result  for  this  same  diode  in  order 
to  give  even  more  insight  into  the  overall  physics  of  its  operation.  Such  further  studies  would  be  valu¬ 
able  for  specific  diode  designers. 


39 


■i 


Appendix  C 


SAMPLE  SOLUTION  FOR  A  VACUUM  TRANSMISSION  LINE 

As  part  of  the  NRL  Light  Ion  Fusion  Research  Program,  high  power  diodes  were  tested  for  use  as 
intense  ion  beam  sources  on  the  AURORA  pulsed  power  generator  at  the  Harry  Diamond  Laboratory. 
For  these  tests,  the  machine  was  operated  in  positive  polarity  so  that  the  magnetically  insulated  coaxial 
feed  line  which  transmits  the  power  from  the  generator  to  the  diode  has  a  configuration  qualitatively 
similar  to  that  pictured  in  the  lower  half  of  Fig.  2.  A  significant  difference,  however,  is  that  in  the 
actual  transmission  line  the  outer  conductor  is  conical  rather  than  cylindrical.  The  specific  geometry  is 
drawn  to  scale  in  Fig.  8  with  all  dimensions  given  in  centimeters.  In  recent  ion  diode  runs  on  this  dev¬ 
ice,  clear  evidence  was  found  that  significant  amounts  of  electron  current  were  impacting  the  outer  cir¬ 
cumference  of  the  end  of  the  inner  conductor  (anode).  This  could  represent  a  significant  loss  of  power 
that  could  have  been  usefully  applied  to  the  generation  of  ions  in  the  diode  instead.  The  2S0  kiloam- 
peres  of  current  flowing  through  the  anode  stalk  had  been  assumed  to  be  quite  sufficient  for  generating 
a  strong  enough  azimuthal  B-field  to  insulate  any  electron  flow  in  the  line.  Apparently  it  may  not  be. 
The  RADBER  code  can  now  be  used  to  check  the  physics  of  the  assumption. 


Fig.  8  —  Modified  AURORA  magnetically  insulated  transmission  line  (MITL) 
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Quite  differently  from  the  case  treated  in  Appendix  B,  the  solutions  at  p  —  ps  proved  to  be 
extremely  erratic  in  the  immediate  vicinity  of  the  true  solution.  This  is  to  be  expected  judging  from 
the  sharp  discontinuity  in  the  slopes  of  both  the  electric  field  and  the  magnetic  field  for  this  case.  On 
the  other  hand,  the  sharpness  of  the  jumps  in  j  space  are  a  useful  aide  in  pinpointing  the  solution  after, 
coarse  tuning  of  the  parameter  guesses  is  accomplished  through  matching  at  p  -  p A.  For  the  3  MV 
case,  the  correct  steady  state  values  for  rs ,  j,  and  (fl«)CJ,h  are  70.51  cm,  3.90  amp/cm2,  and  0.689  kG 
respectively.  Similarly,  for<$0-  4  MV,  they  are  70.33  cm,  4.735  amp/cm2,  and  0.677  kG.  Finally,  for 
<f>o“  5  MV,  the  corresponding  numbers  are  70.03  cm,  6.315  amp/cm2,  and  0.663  kG.  In  Fig.  9  are 
plotted  the  profiles  of  <t>  and  B„  near  the  cathode  for  all  three  cases.  Notice  the  pronounced  discon¬ 
tinuity  in  the  slope  of  each  B„  at  its  respective  electron  sheath  boundary  compared  to  the  relatively 
smooth  transition  of  <t>  across  the  same  boundaries.  In  vacuum,  the  B„  has  a  simple  r~  dependence 
while  <Hr)  is  proportional  to  In  r. 
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The  weakest  insulation  can  probably  be  expected  near  the  end  of  the  anode  sulk  where  the  dam¬ 
age  was  experimentally  observed.  Cylindrically  coaxial  equilibrium  solutions  will  be  sought  for  that 
region  using  RADBER  for  total  voltages,  0O.  of  3.0,  4.0  and  5.0  mV.  Assuming  250  kA  of  axial 
current  in  the  anode  stalk  demands  a  vacuum  Bt  of  approximately  0.7032  kG  at  the  lower  cathode 
radius  of  71.1  cm.  That  value  will  be  used  as  the  scale  B-field,  B0.  RADBER’s  numerical  variables 
must  then  be  set  as  follows: 

a.  LAZ  -  TRUE. 

b.  LC1N  -  FALSE. 

c.  LIONS  -  FALSE. 

d.  ERROR  -  0.01 

e.  RC  -  71.1 

f.  RANODE  -  26.7 

g.  U0  —  3.0E6,  4.0E6,  and  5.0E6 

h.  B0  -  0.7032E3 

i.  CURCEN  -  2.50E5. 

This  problem  is  somewhat  simpler  than  that  of  the  radial  diode  since  no  ion  current  is  present.  Thus, 
only  two  parameters,  j  and  (Bt)alh,  must  be  guessed.  Similarly,  only  two  resultant  parameters,  U  and 
W'  —  Y,  must  be  matched  at  p  —  pj  and  p  —  pA.  As  a  starting  point  for  parameter  guesses,  simple 
Child-Langmuir  theory  predicts  about  6.2,  9.5  and  13.3  amps/cmJ  of  emitted  electron  current  density 
for  0O  -  3.0,  4.0  and  5.0  MV  respectively.  With  such  low  values  for  j,  relatively  little  departure  from 
the  vacuum  Bt  of  0.70  kG  can  be  expected.  This  in  turn  would  lead  to  values  of  ps  very  close  to  pr. 
Because  of  that,  finer  resolution  must  be  used  in  the  numerical  integration  so  that  a  value  of  NDR  — 
750  is  chosen.  The  details  of  the  numerical  solution  search  process  paralleling  that  used  in  Appendix  B 
will  not  be  given  here.  It  is  sufficient  to  note  that  the  reduction  of  the  work  from  a  three-point  to  two- 
point  shooting  problem  permitted  convergence  to  a  solution  for  each  of  the  three  voltages  within  three 
iterations  of  guesses. 
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These  solutions  indicate  that  there  is  very  little  to  fear  from  power  losses  to  electron  flow  in  this 
feed  line.  The  electron  sheathes  at  the  cathode  surface  are  all  on  the  order  of  only  one  centimeter 
thick  compared  to  a  radial  anode-cathode  separation  of  almost  forty-five  centimeters.  The  source  of  the 
experimentally-observed  damage  patterns  is  now  believed  to  be  some  phenomenon  associated  with  the 
plasma  erosion  switches  used  near  the  diode  assembly. 
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